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Abstract. We study the one-dimensional stationary solutions of an integro- 
differential equation derived by Giacomin and Lebowitz from Kawasaki dy- 
namics in Ising systems with Kac potentials, We construct stationary 
solutions with non zero current and prove the validity of the Fourier law in the 
thermodynamic limit showing that below the critical temperature the limit 
equilibrium profile has a discontinuity (which defines the position of the inter- 
face) and satisfies a stationary free boundary Stefan problem. Under-cooling 
and over-heating effects are also studied. We show that if metastable values are 
imposed at the boundaries then the mesoscopic stationary profile is no longer 
monotone and therefore the Fourier law is not satisfied. It regains however its 
validity in the thermodynamic limit where the limit profile is again monotone 
away from the interface. 



1. Introduction 

When hydrodynamic or thermodynamic limits are performed in systems which 
are in the phase transitions regime we may observe perfectly smooth profiles de- 
velop singularities with the appearance of sharp interfaces. We shall study the 
phenomenon in stationary non equilibrium states which carry non zero steady cur- 
rents, the general context is the one where the Fourier law applies, but here it is 
complemented by a free boundary problem due to the presence of interfaces. We 
work at the mesoscopic level considering a model which has been derived in [5] from 
Ising systems with Kac potentials and Kawasaki dynamics, and derive in the hy- 
drodynamic limit macroscopic profiles with an interface which satisfy a stationary 
Stefan problem and obey the Fourier law. 

The mesoscopic model is defined in terms of a free energy functional, the Lebowitz 
and Penrose [L-P] functional (see (|2.1[) in the next section) which is a non local ver- 
sion of the scalar Ginzburg-Landau (or Allen-Cahn or Cahn-Hilliard) functional. 
Its thermodynamic free energy density is obtained by minimizing the L-P functional 
over profiles with fixed total magnetization density and then taking the thermody- 
namic limit where the spatial size of the system diverges. It is found that the phase 
diagram (of free energy density versus magnetization density) obtained in this way 
has a non trivial fiat interval [—mp,mp] (indicative of a phase transition) when 
the inverse temperature j3 is above the critical value (equal to 1 here). This is in 
qualitative and quantitative agreement with the thermodynamics of the underlying 
Ising model with Kac potentials, see Chapter 9 in |5] and references therein. The 
axiomatic theory for such phase diagrams predicts that the values inside (— mp, mp) 
do not appear in any stationary local equilibrium state, so that a macroscopic mag- 
netization density profile will have a discontinuity if it assumes values both smaller 
than —mp and larger than mp. 

l 
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This is just what we see. We fix /3 > 1 and study the stationary solutions of 
dim 

the equations of motion — — = — div /, i.e. div 7 = 0,/ the local current (of the 

conserved order parameter, the magnetization density m here). By a gradient flow 

assumption on its constitutive law, I is supposed proportional to the gradient of 

the functional derivative of the L-P functional: due to the non local structure of 

dfn 

the latter, -— = —div I is an integro-differential equation (see (|2.19[) in the next 

section), which is the same as the one derived by Giacomin and Lebowitz from 
the Ising system, [5], and which has been much studied in the past years, [7], [I], 
[6]. We look for solutions of div I — with a planar symmetry thus reducing to 
a one dimensional problem and prove existence and smoothness of solutions with 
a steady non-zero current. However in the hydrodynamic limit where the size L 
of the system diverges, the stationary profile, once expressed in macroscopic space 
units (i.e. proportional to L), is proved to converge to a discontinuous limit profile, 
solution of a stationary free boundary problem, the stationary Stefan problem, in 
agreement with the axiomatic macroscopic theory. The mesoscopic theory is in this 
respect in complete agreement with the macroscopic one, the mesoscopic profiles 
are smooth versions of the macroscopic ones, they are monotone as well and the 
current is proportional to [minus] the magnetization density gradient in agreement 
with the Fourier law which we may then say to be valid at the mesoscopic level as 
well. 

The mesoscopic theory has however a richer and more complex structure even in 
the macroscopic limit. This is seen for instance if we impose boundary conditions 
which force metastable values at the boundaries, the metastable region being made 
of two separate intervals called the plus and the minus metastable phases (accord- 
ing to the sign of the magnetization) which (together with the spinodal region) are 
contained in the "forbidden region" (—m/3,mp). With boundary conditions one 
in the minus, the other in the plus metastable phases the mesoscopic stationary 
magnetization density profiles are not monotone anymore. We have the "paradox- 
ical" result of a positive [magnetization] current when also the total magnetization 
gradient is positive having fixed at the left and right respectively a negative and a 
positive metastable value of the magnetization. The mesoscopic stationary profile 
is then first decreasing, then increasing and then again decreasing. The Fourier's 
law is therefore not satisfied but, in the thermodynamic limit, the region where the 
profile increases shrinks to a point, which is where the limit profile has a discontinu- 
ity (a sharp interface) . Elsewhere the profile is always decreasing in agreement with 
the Fourier's law (as the current is positive). The stationary profile has therefore 
values all in the metastable region (except at the interface which macroscopically 
is only a point). All the issues presented in this introduction are discussed in some 
more details in the next section, proofs are given in the remaining ones. 



FOURIER LAW AND PHASE TRANSITIONS 



3 



2. Model, backgrounds and main results 

The free energy functional to which we have been referring so far is defined on 
functions m £ L°°(A, [—1, 1]), A a bounded measurable subset of TSL d , as 

Fp t A(m\mAc) = Fp^(m) + - / / J(x,y)[m(x) - m^(y)] 2 

i f >,A(m)= / + \ I I J(x,y)[m(x)-m(y)} 2 (2.1) 

where J(x,y) = J(\x — y\) is a smooth, translational invariant, probability kernel 
of range 1; m^ £ L°°(A C , [—1,1]) is a fixed external profile and 

, , . 1 o 1* \ ci \ 1 + m, /l + m\ 1 — m, /l-m\ 
(pp{m) = --m --S(m), -S(m) = — — log ^— — J + — ^— log ^— J 

(2.2) 

To simplify the analysis we suppose A a cube and consider Neumann boundary 
conditions, namely the functional 

*?TH = jf fl (m) + J | J J nmm (x, y)[m(x) - m{y)] 2 (2.3) 
where J neum (x, y) = J(x, z) with JZ A (y) the set image of y under reflections 

zeR A (y) 

of the cube A around its faces. In d = 1, if A = (r 1 ^!, 1], J noum (a:, y) = J(x,y) + 
J(x, 2(T X — y) + J(x, —leT 1 — y) (e > is a scaling parameter which will vanish 
in the thermodynamic limit). With minor modification what follows in the next 
item "Equilibrium thermodynamics" holds as well for general boundary conditions 
as those considered in (12.11) . 



Equilibrium thermodynamics of the mesoscopic model. 

(The statements in this paragraph are proved in Section 6.1 of [3]). The thermo- 
dynamic free energy density ap(s), s £ [—1,1], is defined as 

ap (s) := \im d inf ^F^ m (m) j m = s j (2.4) 

The limit on the r.h.s. indeed exists and it is equal to: 

ap = cf>*p = convex envelope of 4>p(-) (2-5) 

4>*p = <pp when j3 < 1 and <j>*p ^ <pp when f3 > 1. More precisely let mp be the 
positive solution of 

mp = tanhj^ro/j}, /3 > 1 (2.6) 

then </>|(s), s £ (—mp,mp), is constant and strictly smaller than 4>p(s), while 
4>p(s) — 4>p(s) elsewhere. The values of the magnetization in the interval (—mp, mp) 
are "forbidden" . This is best seen working in the grand canonical ensemble (in other 
words, using Lagrange multipliers). To this end we add a constant magnetic field 
h so that the free energy functional becomes 



Fffi$(m) = F%T{m) -h m (2.7) 



4 



A. DE MASI, E. PRESUTTI, AND D. TSAGKAROGIANNIS 



The grand canonical thermodynamic pressure pp{h) is defined by a minimization 
problem without constraints: 

pp(h) = lim sup f - Fffi%(m) m e L°°(A, [—1, (2.8) 

A— s-B d I M ' J 

Existence ol the limit is again a fact and the two thermodynamics defined by the 
free energy ap and by the pressure pp are equivalent, a property called in statistical 
mechanics "equivalence of ensembles" . Namely pp and ap are inter-related as in 
thermodynamics being one the Legendre transform of the other: 

Pp(h) = sup {/is — ap(s) | s G [— 1, 1]}, ap(s) = sup {hs - pp(h) l/iei} (2.9) 

For any j3 > 1 and any del any maximizcr of ()2.8j) at least for A large enough 
is a constant function equal to vnp^ where mp t h is the solution of the mean field 
equation 

mp. h = tanh{/3(m^ + h)} (2.10) 

which minimizes 4>p{s) — hs and therefore it is not in (—mp,mp), the values in 
(—mp,mp) "are therefore forbidden". 



Gibbsian equilibrium thermodynamics. 

The thermodynamics obtained above is in qualitative and quantitative agreement 
with the thermodynamics of the underlying microscopic model, i.e. the Ising system 
with Kac potential. The Gibbs canonical equilibrium free energy fp.-yim) is defined 
as 

fp, 7 (m) := hm A lun d -^logZ A „ A7 (2.11) 

Z An ,M = J2 X (l E - m )l ^ <5|A„|)e-^- A "^ A ") 

o-a„ e{-i,i} A ™ ieA„ 

where A„ is a sequence of increasing cubes and 

# 7 ,a(cta) = E J 7( a; ^) 'A(a;)crA(y), J 7 (ar, y) = 7 d J 7 (7|x - y|) (2.12) 

(Same free energy is obtained for more general regions and boundary conditions). 
As discussed in Chapter 9 of |9 in d > 2 for any f3 > 1 and 7 > small enough, 
fp,-y(m) is flat in an interval [—mp n ,mp n } and mp n — > mp as 7 — > 0. The 
original result has been proved in [3] and [2] while the fact that in any d > 1, 
lim fp,j(m) = ap(m) is much older and proved by Lebowitz and Penrose, [8]. 



Axiomatic non equilibrium macroscopic theory. 
The basic postulates are (i)-(iv). 

(i) local equilibrium and barometric formula. The free energy of a macroscopic pro- 
file to in the macroscopic (bounded) region fi C ffi d is given by the local functional: 

F™(to) := / ap(m), to e L°°(r!, [-1, 1]) (2.13) 
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(ii) gradient dynamics. The evolution equation in the interior of fi is the conser- 
vation law (D below denoting functional derivative) 

^ = -Vj, j = -xVDFfiT = "XVa^ a'p(s) := ^ (2.14) 

(iii) mobility coefficient, x is a mobility coefficient which depends on the dynamical 
characteristics of the system, we take 

X (s) = 0(1 - s 2 ) (2.15) 

as this is what found when deriving (|2.14p from the Ising spins, [5]-[7]- 

The usual setup for Fourier law has a parallelepiped with different values of 
the order parameter imposed on its right and left faces and Neumann (or periodic) 
conditions on the other ones. By the planar symmetry the problem becomes one 
dimensional and from now on we shall restrict to d = 1 taking 51 = [—£,£]. The 
stationary profiles m(x), x £ (—£,£), verify 

15,3— — = — j — constant, Dp(m) — x(to)o^(to) (2-16) 

and are determined for instance by Dirichlet boundary conditions at ±£, namely 
m(x) — > m± as x — ¥ ±£. To have an increasing profile we shall suppose that 
— 1 < to_ < —171/3 and 1 > m + > mp, the opposite case being recovered by 
symmetry. When j3 < 1 the above is well posed as a'p > but if j3 > 1 the 
denominator vanishes. The macroscopic theory then needs a further postulate: 

(iv) The stationary Stefan problem. There are xq € {—£, £) and j < so that there 
is a solution m(x) of (|2.16[) in (—£, xq) with boundary values to_ and — mp and 

in (xq,£) with boundary values mp and m + . The current — x(m(a;))— a'p(m(x)) 
being equal to j is constant through the interface xq : conservation of mass would 
otherwise impose a motion of the interface against the assumption that the profile 
is stationary. Observe also that since a'a{m) is an increasing function of m in 
[— mp,mp] c and since h(x) is increasing (if j < 0) then m{x) is also increasing. 

A different formulation of the problem is however more convenient for our pur- 
poses. We start by a change of variables, going from m to h. There is a one to 
one correspondence between m and h when {m > mp} and {h > 0} and also when 
{rn < — mp} and {h < 0}. The correspondence is given in one direction by (|2.10l) . 
and in the other by h — a'g(m). Expressed in terms of the magnetic field, (I2.16[) 
becomes 

Hx)= f -^T, m=(a'p)- 1 (h) (2.17) 
Jx a X{m) 

namely m is regarded as a function of h obtained by inverting h — o^(to) and 
X(to) = x(™(^)) becomes a function of h as well. (|2.17|) is then an integral equation 
in h(-) where however xq and j are also unknown: they must be determined by 
imposing the boundary conditions h(±£) = h± := a'p(m±). All this suggests a new 
formulation (alternative to the Dirichlet problem) where we assign Xq and j instead 
of m±. In this way the Stefan problem is written in a compact way as in (|2.17p 
above which is now a "pure" integral equation for h(-) with Xq and j known data. 
We shall mostly use in the sequel this latter formulation when proving that the 
Stefan problem with assigned xq and j can be derived from the mesoscopic theory. 
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As a difference with the Dirichlet problem, in the ll xo,j problem" there is no 
"global existence theorem" , in the sense that given xq and j there are no solutions 
if £ is too large. Indeed (|2.17l) with xq — and j < has a "maximal solution" 
(hj(x),m,j(x)). Namely there is a bounded interval such that 

lirn rrij(x) = ±1, = (2-18) 

(|2.17|) has no solution if £ > £j while any other solution of f|2 . 1 7(1 with the same j 
is obtained, modulo translations, by restricting the maximal solution to a suitable 
interval contained in (—£j,£j). The value lj is strictly finite because the solution 
m(h) of m = tanh{/3/i + /3m} when h — > oo and m(h) — > 1 is to first order given by 
dm 2 dm . . _ 

— -— « p(l -m , Thus — — ~ —j m fl2.16p when m R3 1 hence to(-) converges to 

1 linearly with slope — j (recall j < 0). The collection of all the maximal solutions 
(hj{x),mj(x)) when j G K \ {0} determines in the sense explained above all the 
possible solutions of (|2.17[) . Since lj — > as j — ?• oo and €j — » oo as j — ^ it 
then follows that for any £ the Dirichlet problem with data m± at (m+ 7^ m_ , 
m± in the complement of [— mp, mp]) can be obtained as described above from the 
collection of all the maximal solutions. By taking limits we can also include mp 
and —TO/3. 

By restricting to intervals strictly contained in the maximal interval [—£j,£j] the 
solution (h,m) of (|2.1 7[) is smooth, ||TO,|j < 1, x(to) bounded away from and 
\\h\\ < 00. These are the properties of the macroscopic solution which will be 
repeatedly used in the sequel. 



Stationary mesoscopic profiles. 

Dynamics is defined using the same postulate of the macroscopic theory, namely it 
is the gradient flow of the free energy functional which, in the mesoscopic theory 
is (|2.3p (supposing again Neumann conditions). The gradient flow is (D below 
denoting functional derivative) 

^ = -V7, J = -xV(d^,a) (2-19) 

1 = ~ xv Op iog S " / jncum(x ' v)m[y) dy ) 

With the choice x = /3(1 — to 2 ) (that we adopt hereafter) ()2.19j) becomes the one 
found in [5] from the Ising spins. We suppose again a planar symmetry to reduce 
to one dimension, take A = e _1 [— £,£] interpreting e" 1 as the ratio of macroscopic 
and mesoscopic lengths so that (|2.19[) becomes 

dm = dfdm + m2 d jncum # x 
dt dx V dx dx I 

As in the macroscopic theory it is now convenient to change variables. Define h{x) 

as 

I 1 + TO 

h:=^og-^-J ncnm *m (2.21) 

zp 1 — TO 

Then the current / in (|2.19j) has the expression 

I=- X (m)--, to = tanh{/3 J noum * to + ph\ (2.22) 
dx 
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The stationary problem in the xo,j formulation is then the following. Given any 
xo € (— £, () and j < 0, find m and h so that 

m = tanh{^J noum *m + ^}, h(x) = — (2.23) 

Je-ix X(m) 

We first consider the simpler antisymmetric case where m and h are both odd 
functions. 

Theorem 2.1. Let j ^ 0, xq = 0, £ > and smaller than £j (see (|2 . 18[) ). Then for 
any e > small enough there is an antisymmetric pair (h e (x), m e (x)) which solves 
(|2.23l) in e~ 1 (—£,£) and (h e (e~ 1 x),m e (e~ 1 x)) converges in sup-norm as e — > to 
the pair (h(x),m(x)) solution of the Stefan problem (|2.17[) . Moreover h e and m e 
are both strictly increasing if j < and strictly decreasing if j > 0. 



Remarks, (a) Theorem 12.11 is proved in Section [3] and in Appendix [3] [B] and [C] 
The proof is based on finding the fixed point of the following map: given a function 
h solve the first one in (|2 .23[) to get m and use the second one to find the new h. 
Existence of a fixed point is proved by showing convergence of the iterates h n and 
of the corresponding m n . Since xq = if we start with an antisymmetric function, 
the whole orbit remains antisymmetric and indeed the limit macroscopic solution 
is antisymmetric as well. As we shall see restricting to the space of odd functions 
greatly simplifies the problem. We start the iteration from a profile mo which is 
almost a fixed point: mo is in fact the [scaled by e _1 ] macroscopic solution away 
from while it is equal to the "instanton" (see Section |3]) in a neighborhood of 0. 
We shall prove that all the profiles m n obtained by iterating (|2.23j) are contained in 
a small neighborhood of mo and that the iterates converge to a limit profile m; also 
the corresponding magnetic fields h n are proved to converge to a limit h and the pair 
(h, m) is the desired fixed point which solves (|2.23|) . The crucial point in the analysis 
is to control the change 5m of m in the first equality in (I2.23|) when we slightly 
vary h by Sh. To linear order 5m and 5h are related by {A} ltTn — l)5m — —ph, m 5h 
where Ah t m = Ph,mJ*i the convolution operator with kernel J, and 



cosli {/3J ncum * m + f3h} 
Ph,m = x(m) if m = tanh{/3J ncum * m + f3h} (2.25) 

(the equality Ph.m = x{ m ) m (I2.25[) will be often exploited in the sequel). Thus 
5m = m (—ph,m5h) provided Lh, m := Ah, m — 1 is invertible. In 4 it is shown 
that the largest eigenvalue of Lh, m converges to as e — > and that there is a 
spectral gap bounded away from uniformly in e. By restricting to odd functions 
the leading eigenvalue disappears and the invertibility problem can then be solved. 
As clear from this outline the proof does not give uniqueness which is left open. 

(b) The choice of Neumann conditions simplifies the analysis but other condi- 
tions (provided they preserve antisymmetry) may be treated as well unless they 
contrast with the macroscopic value of the magnetization imposed by j, in which 
case boundary layers may appear, which are instead absent with Neumann condi- 
tions. 

(c) With Neumann conditions the non local convolution term is completely 
defined, but since the evolution involves also derivatives other conditions are needed 
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to determine the solution: our choice was to fix j and xo ■ Dirichlet conditions would 
instead prescribe the limits m± of m[x) as x — > ±e~ 1 £. There are here two types of 
boundary conditions, those which fix m outside the domain and are used to define 
the convolution (in our case replaced by Neumann conditions) and those which 
prescribe the values of m when going to the boundary from the interior (in our case 
are replaced by j and xq). The distinction is not as clear in other models as for 
instance in the Cahn-Hilliard equation where more parameters are involved, we are 
indebted to N. Alikakos and G. Fusco for many enlightening discussions on such 
issues. 

(d) In this paragraph it is convenient to refer to Dirichlet boundary conditions. 
It follows immediately from (|2.14j) that the critical points of the functional are 
stationary, i.e. such that the derivative vanishes, DFp^(m) = 0, A = e _1 [— £, £]. 
They are in fact special solutions of ()2.23|) : those with j = and hence h = 0, 
thus solutions of the mean field equation m = tanh{/3J * (m + toa c )}, when m\c 
is fixed outside A. In this case the limit values of m when x — ¥ dA from the 
interior cannot be prescribed independently, they are generally different from those 
obtained going to dA from the outside by using m\c . If we want different boundary 
values (from the inside) than those produced by solving DFp : h(m) = 0, we must 
look for solutions with a current and we are back to the problem considered in 
this paper. The solutions of the Dirichlet problem with and without currents are 
qualitatively different. In the former case there is a sensitive dependence on the 
boundary values, even macroscopically away from the boundaries while, when the 
current is zero, we see the familiar exponential relaxation towards the stable phases. 

We have a slightly weaker result when xq ^ as we need in our proofs to 
replace the condition h(e~ 1 Xo) — by an integral one, namely / hu* = 0, 

J-e-H 

where u* (whose dependence on e is not made explicit) is a suitable positive function 
on K., symmetric around e~ 1 xo and which decays exponentially as \x — e _1 xo| — > 
oo uniformly in e (if u* were a delta we would then be back to the condition 
h(e~ 1 xo) = 0). We do not control the exact mesoscopic location of the zeroes of 
the magnetization profile m and of the magnetic field profile h, however they differ 
from e~ 1 xo by quantities which vanish faster than any power of e as e — > 0: 

Theorem 2.2. Let j ^ 0, I e (0,£j) and x ^ in (-£,£). Then for any 

e > small enough there is a pair (h e ,m e ) which solves (|2.23j) in e~ 1 (—£,£). 

h e (x e ) = where x e <E e~ 1 (—£,£) and ex c — > xo (see (|G.26|) in Appendix 17?)) . 

Finally (h e (e~ 1 x),m e (e~ 1 x)) —} (h(x),m(x)) in sup-norm as e 0, (h,m) the 
solution of the Stefan problem (|2 . 1 T[) . 



Theorem l2.2l is proved in Section|4]and in Appendix ID ] IE ! IfI and iGl where we derive 
explicit bounds on the speed of convergence. By Theorem 12.11 we can construct a 
quasi solution (/io,mo) of (|2.23l) with an error which around the interface e _1 xo is 
exponentially small in e _1 (we shall exploit this with the introduction of suitable 
weighted norms), (ho, mo) is then used as the starting point of an iterative scheme 
similar to the one in the proof of Theorem 12.11 from which however it differs signif- 
icantly due to the absence of symmetries. The problem is that we cannot restrict 
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anymore to the space of antisymmetric functions and thus need to check that the 
maximal eigenvalue of the operator L obtained by linearizing the first equation in 
(|2.23l) is non zero. We know however from [3] that it is close to zero and actually 
vanishes as e — > 0. But in our specific case we can be more precise and prove that 
it is negative and bounded away from proportionally to e. Thus we can invert L 
but get a dangerous factor e" 1 in the component along the direction of the max- 
imal eigenvector which spoils the iterative scheme as it is and it thus needs to be 
modified. The idea roughly speaking is to slightly shift from e Xo to make smaller 

the component along the maximal eigenvector (hence the condition / hu* = 

J-e- 1 e 

mentioned before Theorem I2.2() and this is enough to make the iteration work. 
The shifts described above are responsible for the derealization of the zero of the 
magnetization profile which may not coincide with that of the magnetic field. 



The Dirichlet problem. 

By Theorem 12.11 and 12.21 it then follows that there are solutions of the stationary 
mesoscopic equation which converge as e — > to the solution of any Dirichlet 
problem with m_ < — mp and m+ > mp or viceversa. At the mesoscopic level, 
though, the boundary values may differ from the prescribed ones but the difference 
is infinitesimal in e. We omit the proof that the above extends to any choice of m± 
in the complement of (—mp,mp) provided m + ^ m_. We thus have a complete 
theory of the derivation of the Stefan problem from (I2.23P gaining a deeper insight 
on the sense in which the values in (—mp,mp) are forbidden. At the mesoscopic 
level in fact such a restriction is absent and in the approximating profiles (h e , m e ) 
which at each e solve (|2.23p . the values in (—mp,mp) are indeed present in m e . 
However the fraction of space where they are attained becomes negligible as e — > 0, 
they concentrate at the interface which in macroscopic units becomes a point and 
in mesoscopic units are described to leading order by the instanton which converges 
exponentially fast to ±mp. 



Under- cooling and over-heating effects. 
In the forbidden interval (—mp,mp) we distinguish two regions: one called "spin- 
odal" is [— m*,m*], to* = y/l - 1//3, the other, {mp > \m\ > m*}, is called 
metastable and it splits into two disjoint intervals, the plus and minus metastablc 
phases according to the sign of m. In the spinodal region <j>p is concave, see (|2.2|) . 
while in (m* , 1) [as well as in (— 1, —to*)] cf>p is strictly convex. If we could restrict 
to (to*, 1) [or to (— 1, —to*)] ignoring or deleting the complement, then <f>p would be 
convex and it could play the role of a thermodynamically well defined free energy 
giving rise to a new "metastable thermodynamics" , new because in the interval 
(rn*,mp) it differs from the "true" thermodynamic free energy ap. When (if ever) 
is it correct to use the metastable one? The usual answer (as its name suggests) 
is that the time scale should not be too long and the initial state of the system 
entirely in the plus [or in the minus] metastable phase. When the evolution is 
given by (|2.19[) initial states entirely in the plus phase (to*, 1) [or in the minus one, 
(— 1, —to*)] evolve remaining in the plus [minus] phase, so that the other values of 
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the magnetization never enter into play and can be ignored. In particular if m e (x, t) 
solves (|2.19|) with initial datum m e (x, 0) = mo(ex), too <E C°°(M. d ; (to*, 1)), then 

lim mJe~ 1 x, <T 2 t) = mix, t) (2.26) 

solution [with initial datum too] of 

— = div ^D*gradmJ, D* = 1 - /3(1 - to 2 ) (2.27) 

with D* > in (m*,l) [and in (-1,-m*) as well]. By ([12]), D* = which 
confirms the interpretation of <f)p as a free energy once we compare D* with the 
expression for Dp in (12 . 16[) . This is proved in [7] where the analysis extends to the 
spin system with Kac potentials, if the Kac scaling parameter is suitably related to 
e so that the time scale is e~ 2 . On much longer times, which scale exponentially 
in e _1 , large deviations and tunnelling effects enter into play with the metastable 
phase becoming unstable, see pQ. 

All the above deals with initial states entirely in the plus [or in the minus] phase, 
much less is known when they coexist. A first answer is provided in this paper, see 
Theorems 12.11 and 12 . 2\ where however the coexisting plus and minus phases are the 
thermodynamically stable ones. In such cases the whole interval {—mp, mp) shrinks 
in the thermodynamic limit to a point, not distinguishing between metastable and 
spinodal values (thus in agreement with the macroscopic, thermodynamics of the 
model) Our next theorem proves that there are also stationary solutions of (|2.23p 
where the plus and minus metastable phases coexist. 

Theorem 2.3. Let j > then for any positive £ smaller than some £j, there is 
an antisymmetric pair (h e (x) , m £ (x)) which solves the stationary problem (|2.23|) in 
e (—1,1) and such that (h e (e x) , m e (e x)) converges in sup norm as e — > to 
(h(x),m(x)) solution of the "metastable" Stefan problem: 

h*(x)=£^, m = <P'p-\h*) m(-£,£)\{0} (2.28) 

h e is strictly decreasing while, to leading orders in e, m e first decreases then in- 
creases (around the origin) and then again decreases. The interval where it in- 
creases has length I e and el t — > as e — > 0. 

The proof of Theorem 12.31 is completely similar to the proof of Theorem 12.11 and it 
is therefore omitted. We did not check that the result extends to the case Xq ^ 0. 
The coexistence of the plus and minus metastable phases is related to the presence 
of a current which "stabilizes" the profile. If j = the stationary solution would be 
close to an instanton except for boundary layers and in the thermodynamic limit 
would converge to the Wulff shape which in this case is simply mp sign(a;). We 
conjecture that the profiles described in Theorem [23] are "metastable" in the sense 
that an additional noise at fixed current j would make the state tunnel toward the 
solution with same j described in Theorem 12.11 
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3. Proof of Theorem 12.11 

In this section we shall prove Theorem 12.11 which will be a corollary of three 
theorems stated below and proved later in three successive appendices. For nota- 
tional simplicity we suppose j < and, as discussed in Remark (a) after Theorem 
12. 1[ we restrict to odd functions, so that by default in this section all functions 
are antisymmetric. The analysis is based on an iterative scheme which is outlined 
in the next two paragraphs. We shall define a sequence (h n ,m n ) which for each 
n satisfies the equality m n — tanh{/3 J ncum * m n + /3h n } and prove that (h n ,m n ) 
converges as n — > oo in sup norm to a limit (h, to) which is the desired solution of 
(I2~23)l . 

The starting element. 

We define h using (|2.21|) with to set equal to to , to the odd function defined for 
x > as 

m (x) = m(x)l[ 0j £ e] (x) + u(e[x - £ e ])l(£^ e - H }(x) (3.1) 

where: to is the instanton (see the paragraph Instanton: notation and properties in 
Appendix [A"f : £ c = x e + 2no, x e : fh{x e ) = rap — e, no a large integer independent 
of e, its value will be specified in the course of the proof of Lemma fA. II as shown 
in Appendix lAl x, scales as loge -1 . Finally, u(r), re [0,£ — e£ e ], is the solution 
of the macroscopic equation (|2.17|) (which in (|2.17|) is denoted by to) . Since ho is 
obtained from mo by (|2.21|) then 

too = tanh{/3 J ncum * to + /3h } (3.2) 

a property which will be satisfied by all the elements of the sequence (h n ,m n ). 
Moreover, denoting by || • || the sup- norm, 

sup||to || < cjxgj < 1 (3.3) 

e 

because ||to|| < mp and ||u|| < 1 since I < lj, see (|2 . 1 8[) and the paragraph 
Axiomatic non equilibrium macroscopic theory in Section [2j 

The iterative scheme. 

As discussed in Remark (a) after Theorem l2.H the idea is to define a transformation 
h — > T(h) [from antisymmetric into antisymmetric functions] in two steps. We first 
find an antisymmetric function to such that to = tanh{/3 J ncum * to + /3/i} and then 
define for x > 

T(h)(x) = -ej [ x(m(v)y\ m = tanh{/3J ncum * to + j3h} (3.4) 

The definition of T{h) thus rests on the possibility of finding an "auxiliary function" 
to which solves the second equality in (|3.4|) and it is such that x(to) -1 is integrable. 
By construction we already know that the auxiliary function too associated to h$ 
exists and ||too|| < cj^j < 1 uniformly in e. The crucial step will then be to prove 
that if h is "close" to ho then (at least for e small enough) there is a unique to "close" 
to too so that the second equality in (|3.4|) is satisfied, ||to|| < 1 and T(h) is thus 
well defined (we do not have general uniqueness as we are in the phase transition 
regime: we cannot exclude that there are other solutions not close to Too). We shall 
then prove recursively that all images h n = T n (h ) are well defined and close to 
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ho, while the auxiliary functions m n are close to mo; moreover (h n ,m n ) — > (h,m) 
in sup- norm as n —> oo. h will then be a fixed point of T with auxiliary function 
m and Theorem 12. II will be proved. 

Notation. 

Our basic accuracy parameter will be e", a € (0,1). e a defines quantitatively the 
a-priori closeness to ho (the elements hk in the iteration will actually be much closer 
to h , \\h k - h \\ < celoge -1 ): 

\\h-h \\<e a , ll/H := sup |/(x)| (3.5) 

\x\<e~ 1 £ 

being understood that all functions we deal with in this section are odd. While 
the basic accuracy parameter clearly depends on e , a € (0,1) above as well as all 
the constants that we shall write in the sequel, denoted by a, b, c and C with or 
without suffixes, will be independent of e. The existence of the auxiliary function 
to in (|3.4p is established next: 

Theorem 3.1. There are constants cj^gj > 1, aj^gj > 0, Cj^-gj := ~jj^p so that for 
all e small enough the following holds. For any h : \\h — ho\\ < e a there is a unique 
m,h in the ball {to : \\m — mo\\ < Cj^^e"} such that — tanh{/?J ncum * m/j + j3h} 
and for any h! : \\h' — ho\\ < e a 

\m h {x)-m h ,(x)\<c m e- a ^ x -y\\h{y)-h'{y)\, x>0 (3.6) 

Jo 

We postpone to Appendix [Al the proof of Theorem 13.11 and proceed with the proof 
of Theorem 12.11 observing that as a consequence of Theorem 13.11 if \\h — ho\\ < e a 
then T(h) is well defined (for all e small enough) because x(to) in the first of p.4[) 
is bounded away from 0. To prove this it suffices to show that ||to| < 1. By (|3.6[) 
with m/j = m and ro/,' = mo, 

||m-rn || < c'^\\h- h \\ < c'^e a (3.7) 

Then by (|3.3I) ||m|| < c j-^i + c '^M e<1 < 1 f° r e small enough. 

Theorem 3.2. There are constants cj^gj and c jg^ > so that for all e small 
enough the following holds. Let ml and to" be both in the ball {m : \\m — mo\\ < 

pX ■ r-X 

c EU ea }> then denoting by h! = ^ -^L, h" = J -^y, 

\h'(x) - h"(x)\ < c m e\j\ f \m'(y)-m"(y)\, x>0 (3.8) 

Jo 

\\hi - h Q \\ <c m e log e~\ h x = T(h ) (3.9) 



We postpone to Appendix [B] the proof of Theorem 13.21 and observe that since the 
transformation T is well defined in the ball \\h — ho\\ < e a we are in business once we 
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show that any iterate of T is in the ball \\h— ho\\ < e a . We postpone to Appendix 
[Cl the proof of: 

Theorem 3.3. There is a constant c r^m > so that the following holds. Suppose 
there is n such that for all k < n, h k — T k (h ) is well defined, \\hk — h Q \\ < e a and 
\\mk — mo|| < e a , mk the auxiliary function in the definition ofT(hk). Then h n is 
well defined and 

Whk+i-hkWKcBzni-yWhx-hol k<n (3.10) 

It is now easy to prove Theorem 12. II We restrict to e > so small that 

Zc&TmcfTmc^e log e -1 < e a (3.11) 

(with C J37T31 defined in (|3.13j) below) and prove by induction that (h k ,m k ) exists 
for all k and moreover \\h k — ho\\ < e a and \\m k — mo|| < e a . Since the statement is 
obviously true for k = we only need to prove that if it is verified for k < n then 
it holds for n as well. By (|3.10j) and (|3.9j) for all k < n, 

\\h k +i - h k \\ < Cf3^oi(^) fc c j 3^eloge~ 1 (3.12) 

which, by (|3.11l) shows that \\h n — ho\\ < e a (for e small enough). Then by Theorem 
13.11 m„ is well defined and by (|3.6|) for all k < n 

\\m k +i - TOfc|| < c (53fy \\h k+1 - /tfcH, CJ333J =max{l,c'jg^j} (3.13) 

Then, using (|3TT2|) , 

||m]fe + i - m k \\ < cj33o}cj333}(i) fe c l (32}eloge _1 (3.14) 

which by p. lip proves that ||m„ — mo|| < e a . Thus the induction is proved and we 
know that for all fc, (hk,rnk) exists, \\h k — ho\\ < e a and \\m k — mo\\ < e a . 

As a consequence of (|3 . 14|) and p,12j) . (h n ,m n ) — > (h,m) in sup-norm with 
h = T(h), m = tanh{/3J ncum *m + fih], and 

\\h — ho\\ < celoge -1 , \\m — mo\\ < ce log e _1 (3.15) 

Making explicit the dependence on e we write the limit as (h e , m £ ) in agreement with 
the notation in Theorem l2.ll Recalling the definition of (ho, mo), see (|3.ip . we then 
obtain the proof of Theorem 12.11 except for the statement about the monotonicity 
of m e which is proved at the end of Appendix [D] 



4. Outline of the proof of Theorem 12.21 

The macroscopic solution. 

For the sake of definiteness we suppose j < and xq > and for notational 
simplicity that the interval (—1,1) is just the interval (—1, 1). By assumption (— 1, 1) 
is then strictly contained in the interval of length 21 j and center xq, i.e. the maximal 
interval where the macroscopic problem with parameters (j, xq) has solution (see 
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the paragraph Axiomatic non equilibrium macroscopic theory in Section [2]). We 
then write I* = 1 + 2xo so that xo is the middle point of the interval [— 1,^*] and, 
for what said above, £* + 1 < 2£j so that the macroscopic problem has a solution 
(ftmacW,"Jmac(i)), % € {—1,1*) with the following properties: it is a smooth pair 
of functions antisymmetric around xq such that ||m mac || < 1 and ||/i ma c|| < °o (so 
that inf x(™mac) > 0). 

The pairs (h*,m*) and (h e ,m e ). 
By Theorem 12.11 for any e > small enough there is a pair (h*(x), m*(x)), x € 
e~ l [— 1,£*] (dependence on e is not made explicit) which solves (|2 . 23[) and is anti- 
symmetric around e~ 1 x Q . Then there is c ^\\ > so that 

P > Ph*,m* > c j4jl for all e > small enough (4.1) 

ft > Ph,m is t rue m general, see (|2.24j) : instead Ph*, m * > c $~T\ because by (|2.25|) 
Ph*,m* = x( m *) = /^(l — ( m *) 2 ) ancl < 1 uniformly in e. This follows from the 

inequality ||m mac || < 1 because, by Theorem l2.1[ lim e ^o H^mac^) — m*(e~ 1 x)| = 0. 
We next define (/i £ ,m e ): 

m e (x) — m*(x), h e (x) — h*(x) + R e (x), xGe _1 [— 1,1] 

R e (x)= J(x,y)[m*(y)-m*(2e- 1 ~y)]dy (4.2) 



We have added the "correction" R e to have: 

m e = tanh{/3[J noum * m e ] + /3h e } (4.3) 

Lemma 4.1. There are rg^J > 0, c ^-^ > and c'g-^j > so that 

ii dm* .. , . dm*(x) . , , 

II —3— II < c E3i' su p I -£ , I ^ c ea e ( 4 - 4 ) 

\x — e~ 1 Xp I >r |4~4] log £ 1 

j4s a consequence \R e {x)\ < cel € -i_i<a;< € -i . 



Proof. By differentiating the equality to* = tanh{/3 J noum <* * m* + /3/i*} (valid in 
the whole interval e _1 (— 1,£*), J noum <* the kernel with Neumann conditions at its 
endpoints) we get 

dm* J rncum,* Jit 

f-^\\<m—, — iiikii + i&ii)<c 

dx dx dx 

because ||ci/i*/dx|| < ce (as h* solves (|2.23[1 ). hence the first inequality in (|4.4I) . The 
second one is not as easy and it will be proved at the end of Appendix [D] Using 
such inequality in (|4.2j) we readily see that |i? e (a:)| < cel e -i_i< a .< e -i, c = cji^j. □ 

Thus R e is "a small boundary field" and except for the small error R e , x( m e)^f 

i i so that the pair (/i £ ,to € ) is "almost a solution" of the stationary problem 
(which could be interpreted as a true solution of a problem with suitably redefined 
boundary conditions). 



An interpolation scheme. 
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A natural way to obtain a true solution from a quasi solution is via the implicit 
function theorem after writing ()2.23[) as a single equation f(h, m) = on the 
space of pairs of L°° functions. Unfortunately we do not have a good control of the 
derivative of f(h, m) which may in principle vanish. The problem simplifies if we try 
to solve only the first one in (I2.23[) and then use the second one to redefine h, which 
opens the way to an iterative scheme as the one used in Section|3] The crucial step is 
the following: find rh such that m = tanh{/3J ncum *m + /?ft,} knowing h and that h is 
"close" to another field h for which there is rh such that m = tanh{/3J neum *m+/3/i}. 
To solve this problem we interpolate writing h(t) = th + (1 — t)h, t € [0, 1], and 
pretending that for all t there is m(t) such that m(t) = tanh{/3 J nenm *m(t)+j3h(t)}, 
we differentiate and get an equation for dm/dt. Its solution will then allow to obtain 



m as m = 




The main point in this procedure is therefore the analysis of the equation for dm/dt. 
This is (|E.2I) in Appendix [El here we just say that it has the form?/; = (Ah, m — ^)~ 1 <t> 
(i/j the unknown), where A h>m = p h , m J ncum *, Ph,m as in (|2.24[) (andp him = x( m ) = 
/3(1 — m 2 ) because m = t&nh{(3J ncum *m + f3h}), J noum * is the convolution operator 
with kernel J noum . The non linearity of the problem reflects in the fact that (h, m) 
above is actually (h(t),m(t)) which is itself unknown but the whole problem boils 
down to an accurate analysis of the operator Ah.m, in a suitably large set of pairs 
(h, m) (the set A in Appendix [D]) . The same problem has appeared in the proof of 
Theorcm l3.ll where however we had the great simplification to restrict to the space 
of antisymmetric functions. In such a restricted space ||^° TO || < 1 for a suitable 
integer n uniformly in e (see Appendix [Al . (1 — Ah^i)^ 1 is then equal to the 
convergent sum J2A^ m and the bound ()A.13j) holds. In the case considered in 
Theorem 12.21 we do not have symmetries and the invertibility of Lh, m '■— Ah. m — 1 
becomes a serious issue. 

In Appendix [Dl we shall establish fine spectral properties of A) lim for all (/i,m) 
in a set A. We shall prove a Perron-Frobenius theorem for Ah. m regarded as an 
integral operator on L°°(e~ 1 [— 1, 1]) showing that it has a maximal eigenvalue A > 0, 
that its eigenvector u (called the " maximal eigenvector" ) has a definite sign (taken 
positive) and, see Proposition ID. 11 that there are positive constants c, d and a so 
that for all e small enough 

0<A<l-ce, 0<m< c 'e- a ^- x °\ (4.5) 

Actually to leading order in e, A = 1 — CjDjjje, see (|D.9j) . A is separated from the 
rest of the spectrum (spectral gap) as stated in Proposition ID.2I 

Our strategy therefore will be to reduce to pairs (h, m) € A, a task accom- 
plished by showing that we can actually reduce to functions h in the very small 
neighborhood Q of h e defined next. 

The set Q 

Let fej4~7j and agjjj be positive parameters (specified in Appendix IF]) . and for any 

/GL 00 (e- 1 [-l,l]) 

( e afT6}(e~ 1 -x) x > e -i X() 

N(f) := sup E e (x)\f(x)\; E e {x) := , +e -,s ~ (4.6) 
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with a^-^ such that d (^^ (xo + 1) = ag^j(l — xq). Recalling that h e is defined 
in (|4.2|) and denoting by u* G L°°(e _1 [— 1,£*],R+) the " maximal eigenvector" of 
Ah* m* we define Q as 



: N(h-h e )<bgj±, J hu* = 

||^f^||<e sup |fcM|< £2 \ (4.7) 

da; |x-e- 1 a; |<(iog€- 1 ) 2 dx > 



The iterative scheme. 
We shall prove in Proposition IF. II that if h G Q then there is m such that m — 
tanh{/?J ncum *TO + /3/i} and moreover (h,m) £ A, A the nice set with good spectral 
properties mentioned earlier. Thus Ah, m has a maximal eigenvalue A with maximal 
eigenvector u, Ah >m u = Xu. In Corollary IF. 41 we shall prove that u is "very close" 
to the restriction of u* to e [— 1, 1], u* the maximal eigenvector of Ah* t m* relative 
to the problem in e~ 1 [— All this collects the properties needed to define the 
iterative scheme and to prove its convergence. We define recursively h n+ \ := T[h n ), 
n > — 1, := h f , as 



h n+1 u* 

h n+1 = h n+1 - " ~ e 1 _ 1 , h n+1 (x) := -ej / xi^iy))^ 1 (4.8) 



(recalling that x( m n) — Ph n ,m n by (|2.25[) 1 . The definition is well posed once we 
prove that h n £ for n > 0, so that there is a unique m„ such that (h n , m n ) £ A. 
We shall indeed prove in Proposition IG.3I that N(h n +i — h n ) < ceN(h n — /i n _i). 
Here we use in an essential way the subtraction in (|4.8p which subtracts [most 
of] the component along the maximal eigenvector u of Af ln _ litnn _ 1 of the "forcing 
term" p/i„_ llTOn _ 1 (h n — 7i n _i). In this way we shall prove iteratively that h n G Q 
so that there is m n with {h n , m n ) G A] moreover we shall see in Appendix [Gl that 

f h * 

h n — > /i and m n — > m as rt — > 00 with m = tanh{/3 j noum * m + h = h — ^ , 

J M * 

h{x) := —ejf e -i x x( m (y))~ 1 - As a consequence the pair (h,m) satisfies (|2.22p 
with /i(x e ) = where x e is such that: 

xMy))" 1 = J ^ ^ (4.9) 

1 x J u 



The proof of Theorem l2.2l will then be completed by showing at the end of Appendix 
[G]that x e exists and that ex e — > xq as e — > 0, see (IG.26I) . 
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Appendix A. Proof of Theorem 13.11 

Before proving Theorem 13.11 we introduce some notation and definitions which 
will be used throughout the whole sequel. 



An auxiliary dynamics 

To construct and compare solutions of m = tanh{/3 J ncum * m + f3h} for given h, 
we introduce an artificial dynamics. Suppose (Ji(t),m(t)), t € [0,1], are smooth 
functions of t and that for all t 



m{t) = tanh{/3J ncum * m(t) + (3h{t)} (A.l) 
By differentiating (|A.1[) with respect to t we get the identity 

dm . dm dh dm dh , . , 

where p t = Ph(t),m(t), Ph,m as in (|2.24[) . and A t = p t J noum *, J noum * the operator 
on L 00 ^" 1 [-£,£]) with kernel J ncum . 

By a change of perspective we now regard (|A.2[) as an equation for the unknown 
^ with p t and considered as "known terms" . We shall prove in this appendix 
that under suitable assumptions on ft, a solution exists and it is unique. We then 

/•* dm 

"construct" m(t) :— m(0) + / and check that it verifies (| A. 1|) . The important 

Jo ds 

point is that the whole procedure works in the same way even if we ask that (|A.1[) 
holds only at time t = 0, being a by-product of the analysis that it remains valid 
for all t € [0, 1]. In the actual applications m(0) = m is a given, known function 
which solves mo = tanh{/3 J ncum * m + f3h }, h(t) = hit + (1 — t)h with h and h\ 
also known and m(t) the unknown, in particular we are interested in its value m\ 
at time t — 1 when h(l) — h%. (|A.2[) then becomes a non linear evolution equation 
and it will be crucial to prove first that L t is invertible, so that the equation can 
be written in normal form 

dm 1 / dh\ . . 

-dt =L * (A - 3) 

and then that L^ 1 ( — Ptjif) is a Lipschitz function of m . 



The operator Ah^ m 

The whole analysis relies on properties of the spectrum of the operator A^ m — 
Ph, m J ncum * (called A t when (h,m) = (h(t),m(t)) as above). We shall study Ah, m 
in a L°°(e^ 1 [—£,£]) setting and since we want to prove that Ah, m — 1 is invertible 
it is crucial to prove that 1 is not in the spectrum of Ah, m - Regarded as an opera- 
tor on L 2 (e _1 [— £,£],p^ m (x)dx), Ah, m is self-adjoint, it has a maximal eigenvalue 
A/i, m which is positive and the corresponding eigenvector Uh >m , called the maximal 
eigenvector, can and will be chosen as strictly positive, see [4]. Further assumptions 
on h and m will allow to prove that Xh.m < 1 — ce, c > 0, and that the rest of the 
spectrum is strictly below 1 uniformly in e. The bound on Xh, m will not be used in 
this appendix, see the proof of (|A.8[) below. 
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Instanton: notation and properties 
The instanton m is a solution of the local mean field equation fh(x) = tanh{/3J * 
fh(x)}, with the following properties (see Section 8.1 and 8.2 of [9]). fh(x) is 

a strictly increasing, antisymmetric function which converges to ±mp as x — > ±oo, 
more precisely there are c ^xa[ an d a ^Ol both positive so that for all x > 

dfh(x) 



We write 



< mp — m(x) < c joi e a 



p — (3(l — m 2 ), A=pJ*, ml 



dx 

dm 
dx ' 



(A.4) 
(A.5) 



where (/)oo = / fP ■ In 0] and Section 8.3 in [9 it is proved that there are 
a {Al§ - > and C ||X6l so that for any bounded function / 

A n (x,y)f(y)dy < ||/|| c S r"H", /=/- (/m') M m' (A.6) 



We can now turn to the proof of Theorem [3J] and restrict hereafter in this appendix 
to the space of antisymmetric functions. After observing that by (IA.4I) 

x.^cloge" 1 (A.7) 

we complete the definition (13.11) of mo by fixing the integer no, chosen so that 

<e- a ^U\l a^>0 (A.8) 

where ip above is any bounded antisymmetric function, recall that ||/|| denotes the 
sup norm of /. Existence of no follows from (|A.6[) because ml and p Xo are symmetric 
and tp antisymmetric so that {ipfh!)^ = 0. 

Lemma A.l. There is a ^AM > so that for any c, a > and all e small enough 
\\A% m il>\\ <e- a ESJ||^||, if\\m-m \\<ce a ,\\h-h a \\<ce a (A.9) 
for any bounded odd function tp. 

Proof. As we shall see (IA.9I) is a straight consequence of (|A.8|) and of 

11^^ -1|| < c 'e a (A.10) 

Pho,mo 

which follows directly from (|3.3[) and the assumptions on h and m. We distinguish 
"small" and "large" values of Xo in A^° m ip(xo)- 

(i). x S [0,x € + n ]. We write 

/ "0 I \ 

i>M \\{ph . mo {^-i)J nmm {x k ^x k ) Ph ' m 7 fc " 1 \ }^ 1 • • • dx nc 
£J{ Pho,m (Xk-l) 

(A.ll) 



FOURIER LAW AND PHASE TRANSITIONS 



19 



Since J ncum has range 1, \xi\ < x e + 2n for all i = 1, . . . n . Then by (lA.7j) for e 
small enough, J ncum (xi,x i+ i) = J(xi,x i+ i). Moreover p h a m ( xi) = p(x % ) (because 
m Q (x) = m(x), h Q (x) = for \x\ < x e + 2n Q ). Thus by (|A 101) 

A n h ° m yj(x Q )~A n ^(xo)\<c'n e a U\\ 

and using (|A.8I) . for all e small enough 

A^(x )\ < e — + c'n e a M < e- a <Mi\\yj\\ 

(ii). xq € [x e + no, e £]. We then write 

A T,m^) = / ^(Xn )l[{ph, m (xk-l)J nCUm (x k - 1 ,X k )dx 1 ---dx no (A.12) 
k=l 

and since J 116 " 111 has range 1, Xi > x e in (|A.12j) for all i = 1, . . .no. When Xi g 
[x e ,£ c ], Pfto.moC^i) = P(^i) < - m(^e) 2 ) and by the definition of x e , m(x e ) = 
771,3 — £• Hence if 6' is such that (3(1 — m^) < 6' < 1 then for all e small enough, 
Ph ,m (xi) < V < 1. When ^ > £ c , p feo mo = /3(1 - tj 2 ) and since u > m/j, 
Ph , mo (xi) < 0(1 - m%) < b' < 1. Thus by (DCTOj) p^a*) < 6 < 1 |^ ro V(^o)| < 
b no \\ip}\. □ 

By (|A.9j) . Lh : m = A} lyTn — 1 is invertible and 



^ = -5>?U fell < i C ™ (A.13) 



where C ||X73l bounds II^mJI- Moreover: 



Lemma A. 2. There exist a ^el > 0, (which defines the parameter introduced in 

(I3.6[) ), c jA.i4t o^d c jA.isb both larger than max{l,- — IA 13t so j or arl y c 

1 - ajAjOi 

and all e small enough 

\L^(x)\ < Cjxnt / e-^GSDl^l^)!, || m - m || < C e a , \\h - h \\ < ce a 
Jo 

(AAA) 

for any x > 0. Moreover if also ml : \\m' — mo\\ < ce a , </ien 

fe - L hL>\ < ctEmWm - m'\\ (A.15) 



Proof. To prove (|A.14I) we write A% m (x,y) as the kernel of A\ m and have by 
(ICTl) 

OO p „ oo 

Lh^W = - E J A tm(x> VMV) = - J E A 1™(X> VMV) 
71=0 n—n(x,y) 
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where n(x,y) > \y — x\ because Ah. m (x,y) = ph, m {x)J nennl (x,y) is supported by 
\x — y\ < 1. (|A.14|) then follows from (|A.9|) . To prove (|A.15|) we write 

use (|A.13|) and that \\A h , m - A h . m <\\ < c\\m - m'\\. □ 



We shall study (|A.3|) with 

h(t) = th" + (1 - t)h f , ti and h" in the ball || h- h \\ < e a (A.16) 

so that \\h(t) — /loll < e a and || - — || < 2e a . The initial datum m' is chosen so that 

dt 

ml = tanh{/3 J ncum * m' + f3h'} and ||m' — mo|| < c'e a where d := /3c ] A.i4 t - To prove 
existence of solutions of (|A.3|) we need to control the "velocity field" 

V(h,m,h) = -L^ m (p h , m h) (A.17) 

where m, h, h are antisymmetric functions. To this end we specify the "free param- 
eter" c which appears in the previous two lemmas so that c > 3c', c' := l3c \&. ul ■ 



Lemma A. 3. For all e small enough the Cauchy problem in the interval t € [0, 1] 

^ = V(h(t)Mt)^), m(0) = m' (A.18) 

has a unique solution m(t) such that \\m(t) — mo\\ < 3c'e a . Moreover, m(t) = 
tanh{/3 J ncum * m(t) + /3h(t)} for all t € [0, 1]. 



Proof. When ||m — mo|| < 3c'e a , \\h — ho\\ < e a the velocity field V(h,m,h) is 
bounded (by fA7[3]) ) and Lipschitz (by ([ATT5]) ). recall that by (|A"7T6|> ||/i|| < 2e Q . 
We thus have local existence and uniqueness till when ||m — mo|| < 3c'e a . Till this 
time 

ll^ll < flU^^U < < 20c^ 
dt dt 1 — aifA~ql 



(recalling from Lemma IA.2I that c j A.i4 t > maxjl, — IA13t — }). Hence ||m(i) — 

i - « rot 

m (0)|| < 2c'e a which ensures existence till t = 1. Recalling (|A.17[) we get from 
(TQg|) that 

dt' 

m(t) — tanh{/3J ncum * m(t) + /3h(t)} is thus constant and being initially it is at 
all times. □ 



■ (m(t) - tanh{/3J ncum * m(t) + ph{t)}^j = 



By taking h' — ho in (|A. 161) by Lemma [A.3l we conclude that for any h : \\h—ho\\ < 
e a there is m which satisfies m = tanh{/?J ncum * m + /3h} and ||m — mo|| < 3c'e°, 
c' = /3c fon . 
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Finally to prove (|3.6j) we write h(t) = th + (1 — t)h' so that 

m h - m h > = J = - J L i r 1 p t (h - h!) 

and (|3.6|) follows from (|A.14j) . The proof of Theorem 13.11 is complete. 



Appendix B. Proof of Theorem 13.21 

By (|3.3[) there is b < 1 so that for all e small enough ||m|| < b < 1 in the 
ball {to : ||m — m || < cj3jjje a }; Q3.8P then readily follows. To prove (|3.9I) we 
observe that h (x) = for x £ [0,x e + 2n ] because in such interval Too = to and 
TO = tanh{/3 J noum * to}. Thus, if hi = T(h ) by (|A~7| 

|/ii(x) — /io(a;)| < celoge , |x| < £ e := x € + 2n 

Define for x > £ e 

xM x(m ) 

then mo(x) — u(e[x — £ c ]) = tanh{(3u(e[x — £ c ] + (3h(e{x — £ c ])}, hence 

m (x) = tanh{pj ncum *m (x)+P(u(e[x-&})-J ncum *m (x) + h(e{x-Z £ }))} (B.2) 
Since m Q (x) = tanh{/3 J neum * mo(i) + Pho(x}, by (JBjJ) 

|fto(a;)-Me[a!-^])|<oe, and by |BT]) |/io(x) — / < ce, 

Since hi(x) = / — ? r, \hi(x) — h (x)\ < \h (x) - / — -|+ce£ £ hence 

Jo x(™o) Je e X(m ) 



Appendix C. Proof of Theorem 13.31 

By assumption for x > 0, < mk(x) < m n (x) + e a 7 k < n. By By Q3.3P ||m || < 1 
so that for all e small enough, Ph k ,m k is uniformly bounded away from 0. There is 
therefore C < oc (recall the current j is a constant) such that 







\h k+1 (x) - h k (x)\ < Ce I \m k {y) - m k -x{y)\ (C.l) 
By (EH) for any y £ [O.e" 1 ], 

\m k (y) - m k -i(y)\ < c I e - a ^\h k (z) - h k ^(z)\ (C.2) 



where we have dropped the suffixes from c and a. We define ip k +\{x) = \h k +i(e x)— 
h k (e~ 1 x)\, x £ [0,£] and by combining (|C.1|) and (|C.2|) we get 

ipk+i{x)<c' ( dy f e-^^-^M^dz (C.3) 
Jo Jo 
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Define Vk(x) — e^ bx tf)k{x), b > a large constant whose value will be specified 
later. We have: 

V k+ i(x)<c' [ e- b{x - y) dy f e- e ' 1 ^ y -^ +b( - z - y) v k {z)e- 1 dz (C.4) 

JO JO 

For e so small that e _1 a > b we have 

K +1 || < c>j\-^ d y ^Kll < C T7 ^- b \HW (C5) 

4c' 1 e _1 2 

We choose b so that — - = -. Then for all e so small that — = < — 

ab 2 e L a — b a 

\\v k+ i\\ < \\\v k \\ which yields ||^ fc+1 || < e b£ (i) fc ||^||. 



Appendix D. Spectral properties of A htTn 

In this appendix we shall first define a set A by weakening properties of the 
pair (h ei m e ) and then prove spectral properties of m when (h,m) is in a small 
neighborhood of A. 

Instanton: additional notation 

Referring to Appendix [X] for definition and properties of the instanton in, we denote 
i _ ,- / i i \ a. i „r i -i_ . 



'x J* 



by m^Q, xo G (— 1, 1), the translate of m by e £0: 

fh xo {x) = fh(x - e^xo), fh' XQ = Px. ( x ) = Pi 1 ~ (a;) 2 ), A Xo :=p 

(D.l) 

Properties of the pair (h ei m e ) 

• m e = tanh{/3J noum * m e + see (|4T3|l . 

• There are r > and 6 > so that p € (x) < e~ b for all \x — e _1 xo| > r. 

din 

• || - 6 1| < c'gjj (proved in Lemma T4. II) and for any c > there is d > so 
that for all e small enough 

sup \m £ (x) - fh Xn (x)\ < cjp^eloge -1 (D.2) 

\x — e~ 1 xq <c log e~ 1 

because by (|3.15j) . ||m* — m || < celoge" 1 . 

• Since — — - = — — and infp £ > then ||/i £ || < c\, II— r— 1| < ce and, by 

dx Pe(x) dx 

.dhJx) —ej /2 , _-. 
sup |— ^| <cie 2 loge 1 



XTie sei A 

By default all coefficients a, c, C with or without a suffix are meant to be positive 
and independent of e; we shall indicate below by item n the n-th property of (ft. € , m e ) 



FOURIER LAW AND PHASE TRANSITIONS 



23 



as listed in the previous paragraph and introduce the quantities (with b in (|D.3|) 
below the parameter entering in item 2) 

QaU > 1 : e-^ElCi-^ioge- 1 = ^ o:=miii{ ^og^OEg} (D.3) 

J = {x : |ac - e _1 a;o| < 2C(HM !ogH}, I' = {x:\x- e -1 £C | < Cfe) loge -1 } 

(D.4) 

(/' will be used later in Proposition ID.1|) . With such notation we define A as the 
collection of all pairs (h, m) such that m — tanh{/3 J ncum *m+0h} and the following 
three inequalities hold: 

Ph,m(x) = 0(1 - m{xf) < e - 2a E3 , \ x - e - l x Q \ > r^ (D.5) 
II -HI < C(5jfr, sup \m(x) - m xo (x)\ < cjgjje log H (D.6) 

(XX X £l 

W<CG3, ll^||<CE3, sup|^-^^|<c E3 e 2 lo ge -(t'-7) 
cfe xe/ dx Px (x) 

where CjdTJ > 2 and: 

• r (C3J > r an d 2a jp~5t = 6/2, b and r are as in item 2. 

• Qp~6i > 2 <%u and c inu > 2c El ( see item 3 )- 

• C fp~7\ > 2max{c, ci} and c fijj\ > 2c[ (see item 4). 
With the above choice of parameters (h e ,m e ) G A. 

Spectral properties in a neighborhood of A 
We continue the analysis of the spectrum of Ah tm started in Appendix lAl assuming 
that (h, m) is in the S ball of A defined as B$(h,m), Bs(h,m) := {(h' ,m') : 

(h.rn)eA 

\\h-h'\\ < 5,\\m-m'\\ < 6}. 

Proposition D.l. There are positive constants C ffijji , c ^9l , c \ d.iq [ > c to.iit an d 
q jp.nt so that for any e small enough there is S — 5(e) > and for any (h, m) in 
the 5 ball of A 

Ph,m > C |D.8t 

(D.8) 

lAfc^-tl-C^l^cgafeloge- 1 ) 2 , = >0 (D.9) 

U TO J /oo 

Moreover, let Uh, m > be normalized as (u\ m )h,m — 1 ^' * s as * n l|D-4|) . i/ien 
sup \u h , m {x) - m' (x)\ < c' mm e(loge- 1 ) 2 (D.10) 

x£l> 

u h ,m(x) < cnmte- aEa|3; - £ ~ lzo1 (D.ll) 



Proof. We shall first prove with slightly better coefficients the inequalities (ID.8[1 
(jD.lip when (h, m) is in A and then use a continuity argument to extend the 
analysis to a 5 ball of A. We thus fix (h, m) € *A and drop the suffix (h, m) when 
no ambiguity may arise. 
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• Proof of (jLX8l> . We bound \m{x)\ < tanh{/3J ncum * ||m|| + /3\\h\\} and \\h\\ < 
CjdTJj hence Ph,m > 2C jp^ , with 2C fpjj = (3(1 — s 2 ), s the positive solution of 
s = t&nh\ (3s + j3C (pj\ } . (|D.8|) then follows in a S ball of (h,m) if 5 is small enough. 

We shall next prove some rough bounds on A and u which will then be improved 
as required in the proposition. We take here (h, m) in a <5-ball of A with S small 
enough. We are going to use repeatedly variants of the obvious equality: 

(fA h , m g) h , m = (fA h ,, m ,g) h ,, m , = [ /J noum * g (D.12) 



We have the lower bound A > - J ° — a °^' m , A = Ay L m and m' here restricted 

((m' X0 ) 2 )h, m 

to A = e~~ [— 1, 1]. Using (|D. 12[) we can rewrite the numerator as 

(m Xo Am Xo ) h . m = / fh' Xo (x)J DCUm (x,y)m Xo (y) = m Xo (x)J(x,y)m Xo (y)- 

JAxA JRxR 

rh' Xo {xf/p X0 + A = f fh' Xo {x) 2 /p Xo + A' 

J A 



= (KJV+ <W 2 — + A', P = Ph,m 
J A PPxa 

where by (fOj) and (|D~3l) , |A| and |A'| are both bounded by < C e~ a (SS e ~ 1 ^~ x ^ < 
ce 2 . The denominator in the last integral is bounded from below because p = 
Ph,m > QrTsl ; r (ID.8[) has already been proved) and p Xo > (3(1 — m 2 ) (as rh(x) 
converges monotonically to mp as x — > oo). By (|D.6j) and for S small enough 
\p(x) — px (x)\ < 2c'jQ-gje log e _1 when x £ I, while in the complement we bound 
fh' Xg as in (|A.4|) (recalling (|D.3[1 ) and use that \p — p xo \ < (3. In conclusion we get 

A> 1-c^rieloge- 1 (D.13) 

with cjeTi3) dependent on C(pj\ , c'jp^ , afpx\ ■ 

• Proof of dDHU- We use (fTXT3T) and the identity u(x) = X~ n (A n u)(x) to get 
upper bounds on u. With n = 1 we obtain 

IMI < A- 1 !! J\\pV2( J u 2 fl 2 < A- 1 !) J\\(3V2( J M^V* < (m4) 

(having used Cauchy-Schwartz and that ||p|| < (3). By tuning n with the distance 
from e _1 a;o we get, using (|D.5|i . 

u(x) < [1 - qrj^eloge- 1 r"e- 2a <13"|M|, when |a; - e _1 a;o| > n + rjp^i 

(D.15) 

which together with (|D.14|) proves (|D.11|) for (h, to) 6 A 

We shall next prove an upper bound on A = A/ lm , (h, to) in a <5-ball of A., S suitably 
small. We start from the operator A xo — p Xo J* acting on _L°°(IR) and since 1 is 

its maximal eigenvalue (with eigenvector m' ), 1 > - — x ° where we choose 

° (w 2 )oo 

f f 

u = Uh^m on A = e [— 1,1] and u = on A c . Denoting (f)oo = / — — , we then 

JR Px 
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have 

{uA X0 u) oo = / u(x)J{x,y)u{y) = j u(x)J neum (x,y)u(y) + R 

J AxA J AxA 

with R = — / u(x)J(x,y)u(y\), y\ the reflection of y into A through its end- 

JAxA<= 

points. By p. lip \R\ < ce -0 ^! 6- ^ 1- * ) ||u|| 2 , hence writing hereafter (•) = {-)h,m, 

(u 2 )*, > (uA^^^iuA^-ce-m^^-^M 2 
= X(u 2 ) - ce- a ^ e ~ 1{1 - xa) \\u\\ 2 

Thus, by (jrXT4l) . A < ^ + C e~ a EB l^^ 1 "^) . We postpone the proof that 

(u 2 ) 

(u 2 ) 

\-rw--l \<ce 2 + ce log e- 1 (D.16) 

and conclude, recalling (|D. 131) and pending the validity of (|D.16I) . 

1 - c frrTa eloge' 1 < A < 1 + g rpTFt eloge' 1 (D-17) 
with c jp~i7ii dependent on Cjd/tJ, c'j^gj, ajn^j. 

Proof of (|rIT6|) 

Recalling that p > C f^M , P < P, Px > /?(! — Trip) and p Xo < /3, we have 

_j, 2 

-JdjH - /,.2 



<■.,:, si < -, ^ < c jmgi (D.18) 



Then, by (|TlTT|) 



u(y) 
(u 2 )U 2 



< ce -*mn\y-^ 1 xo\ i y eA\i' (D.19) 



hence by (|D.3I) and since a fpj[ = — , 



/ M 2 /Px n < C fD-20l (M 2 )ooe 2 , (w 2 )oo> / M 2 /Ps n > (1 ~ C Jp^ol £ 2 ) (m 2 ) qq 

JA\/' Jl' 

(D.20) 



We also have 



(w 2 )oo (-U 2 )oo (u 2 )oo (u 2 )oo 

(« 2 ) 12,1 

7^ 1 < ce 2 + celoge _1 



In the first inequality above we have used (|D. 19|) . in the second (|D.6I) and in the 
third ([TX201) . ([1X16]) is proved. 

We shall next prove (jD.lOj) that we split in an upper and a lower bound for 
u = Uh, m , we take here (h, m) in a (5-ball of A with 6 small enough 
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Proof of (|D.10P (the upper bound). Let y 6 then, writing below yo = y, 

\ n u(y)= f u(y n ) f[{A X0 (y k -uy k ) F^-^M y! ■■■dy n (D.21) 



We choose again n = Cjtoj log e 1 observing that since yo € I' all y k are in /. We 
bound A"" < (l-cjp^aieloge" 1 )-" < (1 + ncelog e - 1 ) < (1 + c'e[log e" 1 ] 2 ). Since 
all y k are in /, by (|D.6|) and for 6 small enough, 

n^p\<i+c e [io ge -T 

hence (with a new constant c) 

/n 
<Vn) \{{A X0 {y k . 1 ,y k )}}dy l ■ ■ ■ dy n (D.22) 
k=l 

We define u so that u(y n ) = (rh'^u)^ m' Xg (y n ) + u. By (| A.6|) - (|D.3|) and for all 

u(y) < KMl 1 + ce(log e - 1 ) 2 ]<m^ t t > oo + ce 2 |M| (D.23) 
which by (|D.14[) can be rewritten as 

<V) <{[l + ce{\o g e-^] { ^^}rn'{y) + c^ (D.24) 



By Cauchy-Schwartz, 

"{!/) ' J 



which, by (|D.16[) . proves 



< + ce (loge- 1 ) 2 ]<(j/)} + ce 2 (D.25) 



<m' (*) + % e( lo g e-Y (D.26) 



( M 2 )!/2 

• Proof of (|D.10I) (the lower bound) . Proceeding in a similar way we get the lower 
bound: 

^ > {[1 - ce(lo g e-^]^^ }<(,) - ce 2 (D.27) 

To bound the curly bracket from below we multiply both sides of (|D.23[) by p~ Q 
and integrate over /'. By ()D.20p : 

(l-qrj^e 2 )^ 2 )^ < « oW ) 2 w [l+ Ce (loge- 1 ) 2 ] + ce 2 loge- 1 || U || 2 
By JEH}, (1 -qrj^e 2 ) < V a ° /o ° [1 + ^(loge' 1 ) 2 ] + ce 2 loge' 1 , hence 



1 - ce(log e- 1 ) 2 < (W ^ U)2 °° < 1 (D.28) 



-^2 ^ \""ao"Voo 
(u 2 /o>. 

,2\ 



which by jEH} yields > (^)V 2 [l _ ce (l og e^f]m' Xo (y) - ce{\oge^) 2 . 

Using (|DTT6|> we then get 



c 



«M,(s) > m' (a?) - -S^oge^ 1 ) 2 (D.29) 
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• Proof of (|D.9|) . We first suppose (h,m) <G A and use for the first time the 
conditions on dm/dx and dh/dx contained in the definition of A. Writing /' for 
the derivative of / w.r.t. x, we differentiate m(x) = tanh{/?J noum * m(x) + /3h(x)} 
and get ml — pj ncum * ml +ph', hence Lm! — —ph', L = A — 1. We multiply both 
sides by p~ x u and integrate over x. Recalling that L is selfadjoint in the scalar 
product with weight we then have 

(A - l)(um') = -(upti) (D.30) 

By (|D.11[) . \{um!) — J p~ x um!\ < ce 2 , having used that \m'\ is bounded, first in- 
equality in (|D.6p . Since ml = pJ nonm * ml + ph' , using the second inequality in 

(EH, ' 

\m'(x) - p(J ncum Y * m(x)\ < sup \ph'\ < ce, x G /' 

y£l' 

Then, by the second inequality in (|D.G|) . 

\m'(x)-pJ ncum *m'(x)\ = \m'(x)-p(J ncum y *m(x)\ < celoge- 1 , x € I' 
and by (jTlTOl) and (|DTTT|) . 

I(um') - (m'm')co| < celoge -1 (D.31) 
Analogous estimates hold for (uph 1 ) and we get 

lA-Il-CEgjell^S^loge- 1 ) 2 (D.32) 

To conclude the proof of the Proposition we need to extend the previous bounds to 
(h, 77i ) in a (5-ball around (/i,m). By (|D~12|) 

A (w 2 ) 

A - ~ (D ' 33) 

The analogous bound can be proved for A/ A and (|D.9j) follows if (5 is small enough. 
The proof of Proposition [Dj] is complete. □ 



The rest of the spectrum is separated from Xh, m by a spectral gap, see j4]. 

Proposition D.2. There are c mja , a flo^^ > 0, cjd35J aTirf a jp^st > so £fta£ 
/or a/Z e small enough the following holds. For any (hi ,m') € A there is 5 = 5(e) 
so that for all (h,m) in a S-ball around (hi ,ml), for all bounded tp 

11^,^11 <qES e " OEBW I^H. V ; = V'~ ^ U H hjn u h , m (D.34) 
\L-£J>(x)\ < CmM J e- a <°m\*-y\^(y)\ d y (D.35) 



The operator A* and its spectral properties. 

We conclude this appendix with a simple extension of the previous results which will 
allow us to complete the proof of Theorem 12.11 and of Lemma FTTl Let (h*,m*) be 
the solution of the antisymmetric problem in e _1 [— 1,£*], with xq the middle point 
in [-1,1*]. We denote by A* the operator p*J ncum <** acting on L°°(e~ l [—l,i*]) 
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with p* — ph* , m * and kernel J ncnm, *(x, y) (defined with Neumann conditions on 
e _1 [— 1,-f*]). We denote by (•)* the integral over e _1 [— 1,1*] w.r.t. the measure 
(p*)~ 1 dx. We first observe that the pair (h*,m*) satisfies the same properties 
(with same parameters) as the pair (h e , m e ) (recall that ra e is the restriction of m* 
to e _1 [— 1, 1] and that h t is the restriction of h* except for the additive term i? e ). It 
then follows that A* and u* satisfy the same properties as Xh.m and u^.m stated in 
Proposition lD.il (without loss of generality we may suppose with same coefficients). 
Also Proposition [D]2] remains valid, indeed its validity is quite general as discussed 
in Section 8.3 of [S]. 

Conclusion of the proof of Theorem 12.11 In order to keep the notation 
used so far we replace the original interval e _1 [— £, £] in Theorem EH] by the interval 
e _1 [— 1,£*] and denote the solution (h e ,m e ) of Theorem l2.1l bv (h*,m*). Reminding 
that it only remains to prove that m*{x) is an increasing function of x (we are 
supposing j < 0), we shorthand ip — and shall prove that ip{x) is strictly 
positive at all x. We have 

^ = L-\-p*^-)=L-\e 3 ) (D.36) 
where L = A* — 1. The positivity of ifc then follows from 

oo 

L-i(ej) = J2(AT(-v) (D.37) 

n=0 

once we prove that the series converges (as all its elements are positive). Conver- 
gence follows because there are a = a(e) and c = c(e) positive such that for all 

n, 

\\{A*) n \\ < ce- an (D.38) 

which would be easy if this was the L 2 norm as we know that A* is the maximal 
eigenvalue and A* < 1 — ce. 

• Proof of (ID.38[) . With A* and u* the maximal eigenvalue and eigenvector of A*, 
u* normalized, ((u*) 2 )» = 1, we have 

(A*) n i> = (X*) n {u*tp)„u* + {A*) n i), ^ = V-<u»*w* (D.39) 
We have A* < 1 - Ce, C > 0, (by (|PT9|> ). we bound u* using ([Dili) , then by (|D~34|) 
IK^TVII ^ C ( A *)"HV'II + c ra e- Q E311"|^|| (D.40) 
hence (|D~38)) . □ 



Conclusion of the proof of Lemma l4.ll It only remains to prove the second 
inequality in (|4T3|) . With V = ^f, by (|D~36)l . and using the previous notation 

i>=([\*-l]- 1 ej J u*)u* + = (ej) - (ej J u*)u* (D.41) 

I [A* - l] _1 ej| < c by (lD~9l) and by (IdTTI) : 

/ u* < c, sup u*(x) < c EI ije~ a l231 r Ei» lose " 1 
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By choosing a fp^ r fz^ > 1 the first term on the r.h.s. of (|D.41j) is bounded by ce 
when \x — £~ 1 %o\ > rjuj loge -1 . The last term is bounded using (|D.34[) by 

oo 

because \\(j)\\ < ce. Lemma |4~T1 is proved. □ 



Appendix E. An auxiliary dynamics 

We return in this appendix to the analysis of the auxiliary dynamics introduced 
in Appendix [XI We shall study the case where initially (ho, mo) £ A and prove 
a local existence and uniqueness theorem under suitable assumptions on h(t). We 
would like to work in A but A itself is not nice in the L°° topology we are using 
as it involves derivatives. For this reason we introduced the <5-balls of A in the 
previous appendix which will play an important role here as well. Our first result 
is a straight consequence of Proposition ID. II and Proposition ID. 21 and its proof is 
omitted: 

Proposition E.l. There is c > and for any e > small enough there is 
S = 6(e) > not larger than the parameter S in Proposition \D.l\ so that for any 
(ho, mo) € A and any (h',m') and (h",m") in the 5-ball of (ho, mo) 

\\L^ m ,\\ <ce-\ \\L^ m ,-L^ m 4<ce- 2 (\\h'-h^ + \\m'-m''\\) (E.l) 

Proposition E.2. Let S and c be as in Proposition \E.1\ and let C be any positive 
number. Then for any e > small enough there is T € (0, 2^7) so mat the following 
holds. For any (ho, mo) G A, and any h(t), t e [0,T], such that h(0) — ho and 
|| ^|| < C there is m(t), t £ [0,T], such that: 

dm 1 . dh\ . s 

~I[ = L t \~ Pt ~ri7h m(0) = mo, L t = L h(t)iTO ( t ), p t = Ph(t), m (t) ( E - 2 ) 

(L t = Lft( t ), m ( t ), p t = Ph(t),m(t))> H m (') _ m °H - I and TO (') = tanh{^J noum * 
to(-) + j3h(-)~\. Finally m(-) is the unique solution of (|E.2|) in ||to(-) — too| < 4. 

Proof. T is determined by the following three conditions: 

T<A ce-^CTK 5 -, (ce- 1 2/3 2 + /?ce- 2 )r< 1 (E.3) 

The first one ensures that \\h(-) — ho\\ < § (because ||^^|| < C); the second one 
(obtained by bounding £j ( — Ptjjf) via Proposition IE. 1[) will imply that ||m(-) — 
mo 1 1 < |, so that (h(t),m(t)) is always in the (5-ball of (ho, mo) and Proposition lE.il 
can be applied. The third condition will imply that the integral version of (|E.2|) 
gives rise to a contraction. 
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Let X := {m £ C([0,T],i oo (e _1 [-l,l]; [-1,1])) : m(0) = m , \\m(-) - m || < -} and 
for m € X let 

V,(m)(i) = mo + f L-' ( - Ps ^l) (E.4) 

By (|E.1|) and the second inequality in (|E.3|> . ||^(m)(t) - m || < ex" 1 f3Ct < -. Thus 

t/> maps into itself. By (IE.1|) and the third inequality in (|E.3[) f/ 1 is a contraction 
with sup norm in x and t. Therefore there is a fixed point m £ X: m = ip(m) and 
since tp maps X into functions which are differentiable in t with bounded derivative, 
m is a solution of (|E.2|) . By (|E.2j) 

(m(i) - tanh{£ J noum * m(t) + /3/i(*)}) = 

so that m(t) - tanh{/3 J ncum * m(t) + /3h(t)} = mo - tanh{/3 J ncum * to + /3/i } = 0. 
Finally, if to solves ()E.2|) and ||to(-) — mo|| < |, then meZ and tjj(rn) = to. Since 
^ is a contraction to is unique. □ 



Appendix F. Properties of the sets A and Q 

The intervals I and i 7 which appear frequently in the sequel have been defined 
in (fLTIj) . 

Fixing the parameters in the set Q 

The coefficient aj^^gj is a positive number strictly smaller than all the parameters 
in Appendix [D] involved with exponential decay, in particular we require a jj^ < 
min{g| fp-5t (1 — x ), a jp-g^ }. The other parameter b f^j\ is fixed so that: 

Q BTl < 1 and such that { j— \ Qfrrt < — 

(F.l) 

Proposition F.l. For all e > small enough if h £ Q there is to shc/i that 
to = tanh{/3 J ncum * m + f3h} and (h, to) G A. 



Proof. Given h £ Q and t £ [0, 1] we define hit) :— th+ (1 — t)/i e observing that 
(/i(0), to(0)) := (/i e , rn e ) € „4 by the definition of A. Let S* be the sup of all s < 1 
such that there exists m(t), t £ [0,s], which solves (|E.2I) in [0,s] starting from 
to(0) = to £ and such that for all such t, (h(t), m(t)) is in the <5-ball of A with 6 as 
Proposition lD.il We shall prove that S — 1 and that for all t < 1 (h(t),m(t)) £ A 
thus proving the Proposition. 

Since || d ^ t - > | = — h t \\ < b ^ji (because h £ Q) we can apply Proposition lE.2l with 
C = ^JTTJ an d (ho, mo) = (h e ,m e ) £ A. As a consequence there is T = T(e) > 
so that m(t) = tanh{^J noum * m(t) + /3h(t)}, t £ [0,T], and 6 = 5(e) so that 
IK*) ~ M < V 2 , \\m(t) - to £ || < <5/2, £ e [0,T]. Since S is not larger than the 
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parameter 8 of Proposition lD.il (sec Proposition [E7T]) we then conclude that S >T. 
By the definition of S the bounds in Proposition lD.il hold for (h(t),m(t)) at any 
t € [0, S] and it is now just a matter of computations to check that (h(t), m(t)) € A 
for all such t. We start by proving that hit) satisfies the conditions in (|D.7|) . 

\\h(t)\\ < \\h e \\ + \\h - h e \\ < + b m < C^g, (as b m < 1< -^p) 

II d Ht) n . ||*e,| || d(ft - fee) „ C{D3 

XjjJU KJUtXJ VAiJb i-i 

for e small enough. Finally in / (defined in (|D.4I0 

< |^i_Z!Z| + i|^_M| < S e2loge -i + e a < c^e^loge- 

dx p xo dx p xo dx 2 

(for e small enough) so that also the last condition in (ID.7|) is satisfied. 

We shall next prove that m{t) satisfies the conditions required in A. We write 
f(t) := — pt[h — h e ]; At, u(t) for the maximal eigenvalue and eigenvector of A t ; (-)t 
for the integral of the measure p^dx on e _1 [— 1, 1]; f(t) := f(t) — (u(t) f (t)) t u(t) . 
By Q 

^ = L^fit) = AT 1 <«(*)/(*))tti(t) + £r7(t), (^W 2 >* - 1 (F.2) 

We bound \f(t)\ < /3N(h - h e )E e (x)- 1 (using that h € Q and p t < /3), u{t) < 
C {dTH e~ a El k-g" 1 ^ I and get 

|<«(*)/(*)>t| < cAT(ft - ^e—BaiC 1 -^)"- 1 

|/(t)| < iV(/i - /i e )( ce - a El}( 1 -™)^ 1 CEn j e -' 1 (llll 2; -^ 1 ™l + ^-^(F.S) 
By dEMJ |L t -7(*)|(») < / e-Eail a! -*l|/(i)|dl/ and 

[ e-amm^-vlE.iy)- 1 < 2 ^ e(x) " 1 

so that, by (|F.2[) and (|F.3[) and since At < ce -1 

|m(t) - m e \ < N(h - h e )( c'e~ 1 e~ a SH {1 " X0 ^ 1 c mli e- a ^ x -^ lx ^ 



which (recalling that N(h — h e ) < b^-ji) proves that for e small enough, 
ii (+\ ii ^ 3/3cjp35ib|[4^ 

||m(t) — m e || < k := — (I 1 .5) 

a$nM - ag^j 

iV(m(t) - m e ) < (Ve^c^nj + 2/?c E^ W ft _ ^) (F . 6) 

By jEI]), Pt < /S(l - (K| - k) 2 ) < e - 4a E3 + 2/3k in {|z - e" 1 ^! > rjc^j} and 
therefore, by (|F.1|) . m(t) satisfies (|D.5[) for £ e [0, T]. To prove the second condition 
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in (|D.6|) we write for x E I, 

\m(x,t) -m xo (x)\ < \m(x,t) - m e (x)\ + \m e (x) - fh xa (x)\ 

< N(m(t) - m^E^xT 1 + cjp^eloge" 1 

< ce^E^xy 1 + CijD^jeloge -1 < 2c^^t log e" 1 

for e small enough (because E^x)- 1 < e^tm^ 1 ^^- 20 ^ l°se -1 ]), The seC ond 
inequality in (|D.6j) then follows recalling that cjjg-gj > 2c |p - 2l . To prove the first 
inequality in (|D.6|) we take the x-derivative of the equality m(t) — tanh{/3J ncum * 
m(t) + Ph{t)}\ 



^ = \-i(u(t)g(t)) t u(t) + Lt 1 ~g( 



9(t) := {- Pt d[k , ke] h m = 9(t) ~ (u(t)g(t)) t u(t) (F.7) 
ax 

By (|4.7p , ||g(i)|| < /3e and an argument similar to the previous one shows that 
|| d y ) || < ce, so that also the first condition in (ID.6|) is satisfied. 

In conclusion we have proved so far that for all e small enough, (h(t), m(t)) € A 
for all i € [0, S]. Suppose by contradiction that S < 1, write S" = min{l, S* + T}, 
then since (h(S),m(S)) € A by Proposition IE.2I there is m(t), t € [5, £?'], which 
solves (|E.2[) in [S, S'] starting from m(S) and such that for all such t, (h(t),m(t)) 
is in the (5-ball of (h(S),m(S)), hence, a fortiori, in the (5-ball of A with 8 as in 
Proposition lD.il This contradicts the maximality of S hence 5 = 1. □ 



Proposition F.2. There are ajp^gj > 0, r jp^ > ; c, c' and a jpi^ > so i/iai /or 
a/Z e small enough the following holds. Let h G Q and (h,m) G A (existence of m 
follows from Provosition \F . 1\) . then 

sup \m(x) - m t {x)\ < e^ 1 ^^ 1 (F.8) 

\x-e- 1 x \<2rfp^\e- 1 

\\h, m - A e | < ce"^ 6 " 1 , ||u fc , m - Ue|| < c'e-^ 6 " 1 (F.9) 



Proof. (|F.8I) follows from (|F.6|) . in the sequel it is convenient to have ajp^gj small, in 
particular a ^r^ < a mjn - Let A and u be the maximal eigenvalue and eigenvector 
of A :— Ah, m , u > normalized so that (u 2 ) — 1 ((•) := (-)h,m)- and A e , u £ the 
maximal eigenvalue and eigenvector of A t :— Af lct m e with u e > normalized so 
that (Mg)e = 1 ((-) e := (-)h e ,m e )- Since (h e ,m e ) and (h,m) are both in A we can 
use the bounds established in Proposition ID. II and ID. 21 for A and A e . 
We then have 

^>Mji .J-'' 1 ; £»>!_. -^n-- (F . 10) 

the first inequality following from (|D. 12|) . To prove the last one we recall that 
Pe = Pfe'.m* > egg, p = ph,m > Cjd^J. In {x : \x - e~ 1 x \ < 2rjF^e~ 1 } we use 
to get 

sup |1 - — | < ce- a fEM e ~ 1 (F.ll) 

\x-€- 1 x \<2r^-gj(- 1 Pe 
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p 2/3 

In {x : \x — £~ 1 xq\ > 2r jp^ e~ 1 } we bound |1 1 < and u e using 

P> < [Tji ( E3 

(jD.111) . Same argument is used to bound from below — - and the first inequality in 

A 

(|F.9|) follows because A and A e are both close to 1 by ce. 



In order to compute the sup in the second inequality in (IF.9I) we consider 
first \x — e xq\ > r fFlH e ■ I n such a case both u and u e are smaller than 
c (|DTT]J e "'^ £ r fEM hence their difference is bounded by c'e -a tEH c , provided 
a firm < Q tD~5l i r i lF~8l - We next take \x — e~ 1 x? l \ < r g\a e ■ Analogously to (|D.21|) 
and with y Q = x, 

X N u(x)= [ u(y N ) TT Wttb-i, Vh) P( ? k - 1 \ }d Vl ■■■dy N (F.I2) 

We choose N — be^ 1 with b > smaller than rjpjj. Then for all k < N, \yk — 
e" 1 xo| < (^ JfTII + b)^ 1 < 2rjFjje~ 1 so that by (|F. 1 1|) for all e small enough, 
u(x) < \~ N [1 + ce- ^ 6 " 1 ]^ A?u(x). We then write Afu(x) = \?(uu e ) e u e + 
A^u(x) so that in {\x — e _1 a;o| < ^ jFTst e -1 } 

u<[l + ce- a ^Y(i^) N (uu 6 ) e u e + X- N \\u\\ Cmm e- a ^ N ) 
By ([Dl4| 

u< [1 + ce-^^Y ((y)* (uu e ) e u e + \~ N ce~ a ^M N ^ (F.13) 

We bound (uu e ) e < {u 2 ) l J 2 , (u 2 ) e = 1 — (tt 2 |l 1) and use the previous bounds 

P< 

for |1 — -2- 1 so that (uu e ) e < 1 + ce~ ae with a and c suitable positive constants. 
By dRTojl ' 

(^l)JV < e -A'io g {i- C e~'"El}^ 1 } < expjc'e-^ ^^ 1 } < I + c^rVES 6 " 1 
A 

Collecting all these bounds and recalling that A < 1 — ce, we get from (|F.13|) 

+ ce - N ( l °s( 1 - c £)- a mm) (f.14) 

hence the upper bound for u in (jF.9[) . The lower bound is proved similarly. □ 



Recall that (h*,m*) is the solution of the antisymmetric problem in e -1 [— 
with xq the middle point in [— 1,£*] and m e the restriction of m* to e [— 1,1], 
We denote by A* and u* the maximal eigenvalue and eigenvector of Ah», m * and 
by A e and u e those of A e = Ah € ,m e , writing u e also for its extension to e -1 [— 1,^*] 
with u e = outside e _1 [— 1, 1]. We suppose (u 2 ) e = ((u*) 2 )* = 1 with the obvious 
meaning of the symbols. 

Proposition F.3. For all e small enough, 

| A* - A e | < ce - a E3J e " I(1 - Xo) , \\u* -u e \\< c j KTae- a EH» £ ' 1(1 ^" ) (F.15) 
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Proof. Sincere = p* ine^-l, 1], (uf), = (u 2 t ) e = 1 so that A* > / u e J ncum -* * u t 
and, by IjD.lip . 

||(jncum,* _ jncum, £ j # ^ < ^-aQ e~ 1 (l-*o) (pjg) 

where J noum ' e and J neum >* are the kernel with Neumann conditions respectively in 
e-^-l.l] and e^-lj*]. Thus 

A* > J u e J nema ' e *u e - ce -«»Ell 6 ~ 1 ( 1 - a "') > A e - ce- a ^^ 1{1 - xo) (F.17) 



For the reverse inequality we write A e > - — — — ^ , the integral being extended 



J u *jneum,e u * 

to e-^-1,1]. Using <|F7l6|> we replace the kernel J neum > £ with J noum <* and then 
extend the integral to e -1 [— 1,^*] bounding u* via (ID. lip which holds as well for 
u* in the whole e _1 [— (see the paragraph "The operator A* and its spectral 
properties" at the end of Appendix |TJ]) . In this way we derive the first inequality 
in (1FT51) . 

As in the proof of Proposition IF. 21 we bound \u*(x) — u e (x)\ < c'e _a EZIl e 
when |x-e _1 x | > rjjrge" 1 using (|D.11|) (supposing a gog < a^ED^Elj)- Wnen 
|x — e a;o| < r l |F^ e~ 1 we write 

u*{x) = (X*)- N (A*) N u*(x) = {\*)- N A?u*(x) (F.18) 

provided (xq + y jFTsl ) e ~ 1 + < e _1 , which is satisfied if N = ae -1 with a > small 
enough. Hence 

= {^) N {u*u £ ) e u e (x) + A?u* (F.19) 

\u*(x) - (^) N (u*u t ) e u.(x)\ < ce~ a ^M N (F.20) 
By ([F~20|) and since by ([PTTT]) (u*) e < c and |((u*) 2 ) £ - 1| < C e" a tEU^ 1(1 " ;ro) 

|1 - (^L)W( u * Ue )2| < cg-aEHjA + gg-oEfle-'a-xo) (F.21) 
A 

so that the second inequality in (|F.15P follows from the first one. □ 
As a corollary of Proposition IF . 2 1 and Proposition IF.3I we have: 



Corollary F.4. In the same context of Provosition [KM, 
|A* - Xh, m \ < ce -"l53l«- I a-"o) i || u * _ u/i m || < cffusie- ^ 6 "^ 1 -^) (F.22) 



Appendix G. Convergence of the iterative scheme 
By (|4.8p with n = —1 we have for x E e -1 [— 1, 1], 

ho(x) = -ej [ xK(y))- 1 = h*(x) (G.l) 
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because m e — m* on e _1 [— 1, 1] and (h*,m*) is a solution of (|2.23|l in e -1 [— !,£*]■ 

Thus by dig) 

f h*u* 

ho (x) = h*(x)- J - T ^- (G.2) 
where the integrals are extended to e -1 [— 1, 1]. Then, recalling (|4.2p . 

/i (a:) - fte(x) = -iJe(s) - V— (G.3) 



Proposition G.l. For all e small enough ho £ Q and 

N(h - h e ) < CE3J e (G.4) 



Proof. / h*u* = because h* is antisymmetric and u* symmetric around 

Je-l (2xq-1) 

the middle point e #0 of the interval e _1 [— (u* is symmetric because the 
eigenvalue A* is simple and A* symmetric). Since the estimates in Proposition lD.il 
apply to u* as well (see the paragraph The operator A* and its spectral properties 
at the end of Appendix |Pj) . by (|D. 1 1|) and since \\h*\\ < c we get 

1 r e-\2x -l) 

u*h* = / h*u* < ce- a lSU e t 1 " 1 "' (G.5) 



Recalling that c p~n in (|4.1[) is strictly positive uniformly in e, we shall next prove 
that 

i 

u * > -S- (G.6) 

-1 C flp.m 



By ([DTl] u* < c nnn , fCTo) then follows from ([4TT]) : 

1 = <(«*) 2 >* = /^ CEm ^{^} / „' 



Thus, recalling (|4.6I) and that a tp~5l i(l — xo) > n fl~6l ; see the paragraph Fixing the 
parameters in the set Q in Appendix IF1 



with < ajo^j < a fp~5i (l — xp) — a ^gt - By Lemma B~T1 N(R e ) < ce which together 
with (|G.7|) proves (|G.4|) . Before proving that ho 6 Q we notice that by (|D.11|) 



|i? e | < cee- a EH e (G.8) 



a property which will be used in the sequel. We have already proved with (IG.4[) the 

first condition for ho € Q- Then, since - — = — it will suffice to show 

dx dx 

that 

I^I^^W^! (G.9) 
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We have 

=^=/ JfryMM-M^-vKdv, ^ = T~ (G.10) 
ax 7c-i cue 

To bound the term IV'(aj) — VK 3 - + £)l> £ = ^ — x ana - • E ' m [ e — ■"•) e_1 + •"■]> 
in (|G.10|) we use the expression (|D.41I) for By (|D.34|) 

||L-V + E( A *)>H < c'^^ES^ < c"e 3 , = (ej) - (ej / u>* (G.ll) 

if iV = Cloge- 1 with C large enough. We have {A*) n cj> = (A*) n (ej)-(ej J u*){A*) n u* 
By (jDlT|) and since A* € (0, 1) 

(A*) n u*(x) < u*(x) < c m jne- a ^ x - e ~ lx ^ 

so that IE^=o{(^*)™^~ ( A *) n ( e J)}\ < ce3 for x > e- 1 -1 and e small enough. 
(|G.9|) will then follow from 

J2 | / (A*) n (x,y)(ej)dy- [ (AT(x',y)(ej)dy\ < ce 2 (G.12) 

n=0 J 

x and x' in [e _1 — 1, + 1]. To prove (|G.12|) we write £ = a;' — x and 

■P*( g»+0 | 



\J(A*) n (x,y)- J(A*) n (x',y)\ < J A*(x, Xl ) ■■ ■ A*{x n - liXn )\l - J[ ^ 



cne6", 6 < 1 



as all points above are in {x : x — e -1 xo > eT l — (N + 1)} (as n < N) and in such 
a region < cgjj < p* < b < 1 (as m* > mp) and \p*(xi + £) — p*(af»)| < ce, by 
Lemma |4~T1 (IG.12|) is thus proved. □ 



Proposition G.2. There are c jg 13 || and cjqTjjj so i/iai /or aZZ e small enough the 
following holds. Suppose that for n > 1, oof/i /i„ and ft. n _i are in Q , then 

N(m n - m„_i) < cjG33j7V(ft,„ - h n -i) (G.13) 

where mi — tanh{/3 j noum * 77^ 4- j3hi}, i = n — l,n. Moreover 

N(m - m e ) < c jqiH e (G.14) 



Proof. We first prove (|G.13[) where we recall that n > 1. Let t £ [0,1] and 
h(t) = th n + (1 — t)h n _i. Since h n and h n _i are in Q then, by convexity, h(t) £ Q 
and by Proposition IF. II there is m(t) such that (h(t),m(t)) £ A, in particular 
m(t) = tanh{/3J ncum * m{t) + /3h(t)} and m(0) = m n _i, m(l) = m„ so that 

|m„ - m n _i| < sup |. By (|R2l) . (|DTll and (jL\35]) . recalling that p t < /3 

te[o,i] 

and writing 

^(ar) := f | / «(*)[&„ - /i„_i]f) e^EUl*-*" 1 ^! (G . 15) 



^ce-^j + c I e- a ^m^(\h n -h n ^\(y)+^(y))dy (G.16) 



rff 
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We are going to prove that 

1> < cfGJne w e- a ^ x -^ lx ^N(h n - h n -i) (G.17) 



By the definition of Q, J u*hi = 0, i = n — 1, n, then 

l{t)[h n - h n -i] = [ [u{t) - U*][h n - /l n _i] (G.18) 



| J [ U (t) - U*][h n - h n -x}\ < N{h n - hn-x) I \u(t) - U*\E^ < ce 10 N(h n - hn-x) 

(by Corollary IfH]) . (|G.17|) is proved. Using (|G.17j) we have 

j e- a mB\x-y\^y)dy<ce w e- a \ x -^ lx °\N{h n -h n -x), a = nunjaflQl, a^m\ 

The other integral on the r.h.s. of (|G.16[) is bounded by 

J e -*mm\*-y\\h n -h n -x\{y)dy < N(h n -h n -x) j ' e- a ^M^E e {y)dy 

< ce - a EM\ x ~ <L ~ lx °\ 
because ajp-jsj > ag-gj. Collecting all these bounds we have from (|G.16p 

\^^-\{x) < c^e^Ell 1 -' -11 "! + (1 + e 10 )e- a ^ x -^ lxo ^N(h n - h n -x) 
which proves (IG.13I) . 

The proof of (|G.14j) goes in the same way except for (|G.18j) which becomes 

(t)[h -h e }= [ [u(t) - u*][h -h e ]+ [ u*h t (G.19) 



By (|G.5j) and (|G.8|) the latter integral is bounded by < ce -0 ^ 6 ^ 1-2 ^ and the 
bound (|G.14[) is not affected. 

□ 

Proposition G.3. There is c jQ 2 p|| > c Hga[ so that for all e small enough the 
following holds. Given any n > if hk, k < n, is well defined and in Q then also 
h n+ x is well defined and 

N(h k+ x-h k )<l C ^ eN(hk - kk -^ k = 1 - H (G.20) 



Proof. By Proposition IF. II there is m k , < k < n, so that (hk,rrik) & A and, by 
> Cjp~8t ■ As a consequence p„ is bounded and h n+1 is well 
defined; moreover — p^T-il — c \ m k — Ta k -\\ and (for x > 6~ 1 xq) 

( f u *f\ r 

\hk+x - h k \ < eel/ + r , f(x) = \m k -m k -x\dy 

Let x > e^xo, then by (IG.13|) for k > 1 
e a E21 ( e - 1 -,) /(x) = | f" e - a(m (x-y) cm:m}N ( hk _ hk _j < cN r hk _ hk _j 

(G.21) 
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and by ([Pll]) 




< cN(h k - h k -i)e~ a ^ e 1(1_:Eo) (G.22) 

By (|G.6I) J u* is bounded away from hence the bound in (|G.20[) for k > and 
x > e^xo. When A: = we use (|G. 14[) after bounding |too — m e | < N(mo~ m e )E~ 1 . 
Analogous bounds hold for x < e _1 a;o and (|G.20[) is proved. □ 

Proposition G.4. In the same context of Provosition I G. Si for any k < n + 1 

N{m k - m e ) < ce, N(h k - h e ) < c'e (G.23) 

where c = cj^nj + , c = c fQ^ (l + ). Moreover 

1 - ec tO20l 1 - ec jG^Qi 

da; 

and, m particular, h n+ i E G ■ 



Proof. By (10201) for i > 0, iV(^+i - h) < (ec(gm) i+1 and by (lG~4l) . 7V(/i - 
fr-i) < ecjcjj < ec jG^ol , h-i = h e . Then 

iV(fe fc - h e ) < V JV(hi - fej-i) < ecfiT2ot(l + ) 

~£ — 1 ^ ec fsm 

hence the statement in (|G.23|) about hk- The one about nik is proved similarly, 
using flSZHJ) and (|GT3|) . To prove ([GT24)l wc write 

d(h fc - /i fc -0 < ^ _ m < d ^ _ hk ^ 

ax 

so that by (|G.9|) . (|G.20|) and (IG.4|) and with h-t := /i £ , for a; > e _1 x , 
ax ax 

i—l 

An analogous bound holds for x < e~ 1 xo hence (|G.24p . □ 



Conclusion of the proof of Theorem 12.21 We shall first prove by induction 
that h n £ Q for all n. Indeed ho £ Q by Proposition IG . 1 1 and by Proposition IG. 31 if 
hk £ Q for all k < n, then h n +\ £ Q. Thus h n £ Q for all n and by Proposition lF.il 
there is m n so that (h n ,m n ) £ A. We shall next prove that (h n ,m n ) converges in 
sup norm to a limit (h, m) and that, writing h—\ = h e and m_i = m e , 

oo oo 

h = h e + y~](hn - h n -x), m = m e + ^(m„ - m„_i) 

n=Q n=0 

The first series in fact converges because N(h n+ i — h n )< (c j G 2 o t e)" +1 , as remarked 
in the proof of Proposition IG.4I The series for m converges for the same reason 
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because N(m n — m„_i) < c \G.i3\ N(h n — h n -\). By ()G.23j) . N(m — m e ) < ce and 
N (h — h t ) < ce; moreover 

m = lim m„ = lim tanh{^J noum * m n + /3/i n } = tanh{/3J noum *m + /3h} 

n—too n—^oo 

fh n u* - f x -ej 

because h — lim {h n }, h n (x) = / — -. As a consequence, for 

n ^°° J u* ' Je-ixo X{m n ) 

any z £ e _1 (-l, 1), 

nW. 

-ej 
X(m) 

so that the proof of Theorem 12.21 will be complete once we show that: 



h(x) = h(z) + 



• there is x t such that h(x e ) = • lim ex e = x$ 

£-10 

The existence of x e is proved using the implicit function theorem. We thus want 
to prove that h(e~ 1 xo) is "small". Since /i(e _1 Xo) = we need to control | / hu*\. 
We write J hu* = f(h — h n )u* + J (h n — h )u* + J h u*. The first term vanishes 
as n — > oo because J u* < oo and for any x > e~ 1 xg (for instance) 

\h(x)-h„(x)\ < \ej\ IxO)" 1 - x{ m n-iY 1 \ < ce\x - e~ l x \\\m - m„_i|| 

Je~ 1 x 

< c \\m — m n _i|| — > 

having used that x( m n) = Ph n ,m n > C jPTa an d therefore x(m) > Cmjj as m n — » m 
in sup norm. Analogously |x(w„_ 1 ) _1 — x( m e) _1 | < c\ m n-i — w e | < cN(m n _i — 
m^E^ 1 < deEj 1 , so that for x > e _1 a;o 



M»)-'H)(«)I < Itfl / lx(m n - 1 )- 1 -xK)" 1 |<ce 2 / ^(y)" 1 

«/ e 1 a:o ** e 1 xq 

< ce 2 \x-e- 1 x \e- a ^i e " 1{1 - x ' ) 

Thus by (jD.llj) the second term is bounded by J |fr„ - /t |u* < ce 2 e -o EH( 1-a!0 ) e_I . 

Finally, since by (lG~T) fr = h*, by ([G~5]) y h w* < c e- a EU e In 
conclusion, letting n oo, 

We shall next prove that h is continuous and that it changes sign in a small in- 
terval around e xq, thus concluding that there is x e in such interval where h 

vanishes. We have —(a;) = —, \, which, by (|D.8I) . is bounded. Moreover 

ax p(l — m (X)) 

N(m — m c ) < ce and \m e — fh XQ \ < ce in [e -1 £o — 1, e _1 xo + 1]. In such interval 

therefore 1^— (x)\ > ae, a > 0. Hence there is x t where h(x e ) = and 
ax 

\xe - e _1 x | < c"ee- a ^ < - 1 - Xo)e ~ 1 (G.26) 
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